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Magnetic reconnection is a process of magnetic 
field topology change, which is one of fhe mosf 
fundamental processes happening in magnetized 
plasmas. In most astrophysical environments the 
Reynolds numbers corresponding to plasma flows 
are large and therefore fhe transition to turbulence 
is inevitable. This turbulence, which can be pre¬ 
existing or driven by magnetic reconnection itself, 
must be taken into account for any theory of magnetic 
recoimection that attempts to describe the process in 
the aforementioned environments. This necessity is 
obvious as 3D high resolution numerical simulations 
show the transition to the turbulence state of initially 
laminar reconnecting magnetic fields. We discuss 
ideas of how turbulence can modify reconnection 
with the focus on the Lazarian & Vishniac (1999) 
reconnection model. We present numerical evidence 
supporting the model and demonstrate that it is 
closely connected to the experimentally proven 
concept of Richardson dispersion/diffusion as well 
as to more recent advances in understanding of the 
Lagrangian dynamics of magnefized fluids. We point 
out that the Generalized Ohm's Law that accoimts 
for turbulent motion predicts the subdominance of 
the microphysical plasma effects for reconnection for 
a realistically turbulent media. We show that on 
of the most dramatic consequences of turbulence 
is the violation of the generally accepted notion 
of magnetic flux freezing. This notion is a corner 
stone of most theories dealing with magnetized 
plasmas and therefore its change induces fundamental 
shifts in accepted paradigms, for insfance, turbulent 
recoimection entails reconnection diffusion process 
that is essential for understanding star formation. 
We argue, that at sufficiently high Reynolds numbers 
the process of tearing reconnection should transfer 
to turbulent reconnection. We discuss flares that 
are predicted by turbulent reconnection and relate 
this process to solar flares and gamma ray bursts. 
With reference fo experiments, we analyze solar 
observations, in-situ as measurements in the solar 
wind or heliospheric current sheet and show the 
correspondence of data with turbulent reconnection 
predictions. Finally, we discuss First Order Fermi 

arcpleratinn nf nartirles that is a natural rnnspniipure 
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1. Problem of Magnetic Reconnection in Realistically Turbulent 
Plasmas 

Magnetic fields are known to critically modify the dynamics and properties of magnetized 
plasmas. It is generally accepted that magnetic fields embedded in a highly conductive fluid 
retain their topology for all time due to the magnetic fields being frozen-in [1,2]. This concept 
of frozen-in magnetic fields is a basis of many theories, e.g. of the theory of star formation in 
magnetized interstellar medium. 

In spite of this, there is ample evidence that magnetic fields in highly conducting ionized 
astrophysical objects, like stars and galactic disks, show evidence of changes in topology, i.e. 
"magnetic reconnection", on dynamical time scales [3-5]. Historically, magnetic reconnection 
research was motivated by observations of the solar corona [ 6 - 8 ] and this influenced attempts 
to find peculiar conditions conducive for flux conservation violation, e.g. special magnetic field 
configurations or special plasma conditions. For instance [5] showed examples of magnetic 
configurations that produce fast recoimection and much work has been done showing how 
reconnection can be accelerated in plasmas with very small collision rates [9-13] (see also reviews 
[14—16] and references therein). However, it is clear that reconnection is a ubiquitous process 
taking place in various astrophysical environments. For instance, magnetic reconnection can be 
inferred from the existence of large-scale dynamo activity inside stellar interiors [17,18]. We would 
argue that it is also required to enable the eddy-type motions in magnetohydrodynamic (MHD) 
turbulence, e.g. in the Goldreich & Sridhar turbulence [19]. In fact, it is easy to show that without 
fast magnetic reconnection magnetized fluids would behave like Jello or felt, rather than as a 
fluid. 

It is clear that solar flares [20] are just one vivid example of reconnection activity. Other 
dramatic reconnection events attributed to recoimection include 7 -ray bursts (see [ 21 ] for a 
review), while reconnection routinely takes place in essentially everywhere both in collisional and 
collisionless magnetized plasmas. Incidentally, magnetic reconnection occurs rapidly in computer 
simulations due to the high values of resistivity (or numerical resistivity) that are employed at the 
resolutions currently achievable. Therefore, if there are situations where magnetic fields reconnect 
slowly, numerical simulations do not adequately reproduce the realities of astrophysical plasmas. 
This means that if collisionless reconnection is the only way to make reconnection rapid, then 
numerical simulations of many astrophysical processes including those of the interstellar medium 
(ISM), which is collisional, are in error. Fortunately, observations of collisional solar photosphere 
indicate that the reconnection is fast in these environments (see [ 22 ]), which contradicts to the 
idea that being collisionless is the prerequisite for plasma to reconnect fast. 

What makes reconnection challenging to explain is that it is not possible to claim 
that reconnection must always be rapid empirically, as solar flares require periods of flux 
accumulation time, which corresponds to slow reconnection. Thus magnetic reconnection should 
have some sort of trigger, which should not depend on the parameters of the local plasma. In this 
review we argue that the trigger is turbulence. This opens a wide avenue for the application of 
turbulent reconnection theory to explain the astrophysical explosions, e.g. solar and stellar flares 
and superflares, as well as gamma ray bursts. 

A lot of support to models of reconnection based on plasma physics comes from the in situ 
measurements of magnetospheric reconnection. While important for some practical purposes, 
e.g. for some aspects of Space Weather program, this reconnection happens on scales comparable 
to the ion inertial length and therefore is atypical for large scale reconnection that happens in 
most astrophysical systems. We argue that the large scale magnetic reconnection is based on 
MHD turbulence physics making small scale plasma reconnection processes irrelevant for the 
reconnection rates that are attained. 

With the advent of numerical simulations it becames clear that the regular schemes of 
reconnection, like classical Sweet-Parker or Petschek reconnections, do not work. Instead, the 
reconnecting systems transfer to a more chaotic state which is characterized by the hierarchy 
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of magnetic islands in 2D [23-28] or a more complex chaotic state in 3D [29]. We argue that as 
the scale of reconnection layers increases the turbulent reconnection will inevitably take over, 
modifying and suppressing the plasmoid instability that gives rise to currently observed picture. 

Turbulence generation has long been associated with magnetic reconnection processes (see 
[30]). This review, however, is mostly dealing with how turbulence changes the rates of magnetic 
reconnection, although we also consider turbulence generation by reconnection. 

Magnetic reconnection is a ubiquitous process in turbulent media, but it is not easy to observe 
as reconnection transfers most of fhe energy into kinetic motion related to smaller scale eddies, 
thus supporting the energy cascade. Apart from Solar flares, which dynamics can be compared 
with the predictions of turbulent reconnection, in-situ measurement of reconnection available for 
the solar wind provide ways of testing theoretical predictions. We show that both sets of data are 
consistent with turbulent MHD based magnetic reconnection. 

It is worth noting, that our discussion addresses 3D magnetic reconnection. The change of 
dimensionality of physical problems changes frequently the physics involved. For the theory 
based on MHD turbulence, this is very important to note that the properties of MHD turbulence 
are very different in 2D and 3D. 

The theory of turbulent reconnection that we describe is based on the Lazarian & Vishniac 
work ( [31], henceforth LV99) and the extensions of the original model in subsequent publications, 
for insfance in Eyink, Lazarian & Vishniac ( [32], henceforth ELVll). The original LV99 model 
was supported by numerical simulations, some results of which have been published (see 
[32-34]), as well as different pieces of observational evidences that we describe in the review. 
Additional theoretical support for the model comes from very recent work by [35]. While our 
review reflects our optimism based on the successes of the LV99 model in explaining different 
astrophysical phenomena, e.g. gamma ray bursts ( [36] and ref. therein), removal of magnetic 
fields from molecular clouds (e.g. [37] and ref. therein), we feel that the challenges presented 
by the variety of astrophysical conditions are very stimulating for further studies of magnetic 
reconnection. We also accept the limitations of our model that is intended for describing the 
astrophysical phenomena at large scales and therefore adopting MHD approximation. Therefore, 
magnetic reconnection happening at the scale of ion Larmor radius, as is the case of fhe Earth 
magnetosphere cannot be described by the model. Important cases of magnetic reconnection 
in the presence of plasma effects as well as plasmoid instabilities are described in an extensive 
review by [16]. Other cases when magnetic reconnection can be fast in MHD regime without 
turbulence are discussed at length e.g. in an excellent book by [38]. Thus this review should be 
viewed as a personalized outlook on the reconnection problem by the authors who are exploring 
the connection of fhe two ubiquitous processes, namely, magnetic reconnection and astrophysical 
turbulence, while many issues of fhe problem are far from being finally settled and different ideas 
and alternative models are being tested and explored by different research groups. We accept that 
magnetic reconnection, similar to magnetic turbulence, is a very deep subject where the synergy of 
different approaches and techniques may prove to be beneficial eventually. We also note that the 
claim that turbulence can accelerate magnetic reconnection predates the LV99 model (see [39-43]). 
Some new approaches to turbulent reconnection were formulated more recently (see [44]). In the 
review we provide a comparison of LV99 with these approaches. 

In what follows, we argue that turbulent reconnection is the generic process taking place 
in astrophysical environments which are turbulent due to the huge Reynolds numbers of the 
flows involved. The turbulence can be pre-existing and also self-generated by the recoimection 
process. We provide the MHD description of astrophysical turbulence in section 2, describe 
LV99 model of turbulent reconnection in section 3, provide its elaboration and extension in 
section 4, demonstrate examples of the numerical testing of the model in section 5, discuss the 
observational testing of the model with solar data and solar wind data in section 6, outline the 
implications of the model in section 7, provide a comparison of fhe model with other ideas of fast 
stochastic reconnection in section 8. We conclude by discussing the general tendency of models 
of reconnection to get more stochastic in section 9. 
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2. Astrophysical Turbulence and Its MHD Description 


Observations of the interstellar medium reveal a Kolmogorov spectrum of electron density 
fluctuations (see [45,46]) as well as steeper spectral slopes of supersonic velocity fluctuations 
(see [47] for a review). Measurements of the solar wind fluctuations also reveal turbulence power 
spectrum [48]). Ubiquitous non-thermal broadening of spectral lines as well as measures obtained 
by other techniques (see [49]) confirm that turbulence is present everywhere in astrophysical 
environments where we test for its existence. This is not surprising as magnetized astrophysical 
plasmas generally have very large Reynolds numbers due to the large length scales involved and 
the fact that the motions of charged particles in the direction perpendicular to magnetic fields are 
constrained. Laminar plasma flows at these high Reynolds numbers are prey to numerous linear 
and finite-amplitude instabilities, from which turbulent motions readily develop^. 

Indeed, observations show that turbulence is ubiquitous in all astrophysical plasmas. The 
spectrum of electron density fluctuations in Milky Way is presented in Figure 1, but similar 
examples are discussed in [48,53] for solar wind, [54] for molecular clouds and [55] for the 
intracluster medium. The plasma turbulence is sometimes driven by an external energy source, 
such as supernova in the ISM [56,57], merger events and active galactic nuclei outflows in the 
intercluster medium (ICM) [58-60], and baroclinic forcing behind shock waves in interstellar 
clouds. In other cases, the turbulence is spontaneous, with available energy released by a rich 
array of instabilities, such as magneto-rotational instability (MRI) in accretion disks [61], kink 
instability of twisted flux tubes in the solar corona [62,63], etc. In all these cases, turbulence is 
not driven by reconnection. Nevertheless, we would like to mention that an addihonal driving 
of turbulence through the energy release in the reconnection zone can sometimes be important, 
especially in magnetically dominated low beta plasmas. We discuss the case of turbulence driven 
by reconnection in section 4c. All in all, whatever its origin, the signatures of plasma turbulence 
are seen throughout astrophysical media. 

As turbulence is known to change dramatically many processes, in particular, diffusion and 
transport processes, it is natural to pose the question to what extent the theory of astrophysical 
reconnection must take into account the pre-existing turbulent environment. We note that even if 
the plasma flow is initially laminar, kinetic energy release by reconnection due to some plasma 
process, e.g. tearing and related plasmoid generation, is expected to generate vigorous turbulent 
motion in high Reynolds number fluids. 

Turbulence in plasma happens at many scales, from the largest to those below the proton 
Larmor radius. To imderstand at what scales the MHD description is adequate one needs to 
reiterate a few known facts [32,64]. Indeed, to describe magnetized plasma dynamics one should 
deal with three characteristic length-scales: the ion gyroradius pi, the ion mean-free-path length 
£mfp,i arising from Coulomb collisions, and the scale L of large-scale variation of magnetic and 
velocity fields. 

The MHD approximation is definitely applicable to "strongly collisional" plasma with 
£mfp,i ^ Pi- This is the case, e.g. of star interiors and most accretion disk systems. For 
such "strongly collisional" plasmas a standard Chapman-Enskog expansion provides a fluid 
description of the plasma [65], with a two-fluid model for scales between £mfp,i arid the ion 
skin-depth 5i = pi/y/Jk and an MHD description at scales much larger than 5i. 

Hot and rarefied astrophysical plasmas are often "weakly collisional" with £mfp,i Pi- 
Indeed, the relation that follows from the standard formula for the Coulomb collision frequency 
(e.g. see [66]) is 


£mfp,i A V[4 
Pi * InA c ’ 


( 2 . 1 ) 



where A = A-kuX^ is the plasma parameter, or the number of particles within the Debye screening 
sphere. For some media that A can be large. 


MIT addition, the mean free path of particles can also be constrained by the instabilities developed on the collisionless scales 
of plasma {see [50-52]). In this situation not only Alfvenic but also compressible turbulent modes can survive. 
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Figure 1. Big power law in the sky from [45] extended to scale of parsecs using WHAM data. From [46]. 


For the "weakly collisional" but well magnetized plasmas one can invoke the expansion 
over the small ion gyroradius. This results in the "kinetic MHD equations" for lengths much 
larger than pi. The difference between these equations and the MHD ones is that the pressure 
tensor in the momentum equation is anisotropic, with the two components py and px. of 
the pressure parallel and perpendicular to the local magnetic field direction [64]. In "weakly 
collisional", i.e. £.^fpi., and collisionless, i.e. £rnfp,i^ L systems turbulence is bound to 
induce instabilities that limit the effective mean free path [£rnfp,i\eff by magnetically mediated 
scattering of particles [51,67]. This effective mean free path is a game changer and it is not 
surprising that numerical simulations in [68] that accounted for this effect demonstrated that 
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turbulence in "collisionless plasmas" of galaxy clusters is very similar to MHD turbulence on the 
scales larger than [emfp,i]eff- 

We can also note that additional simplifications that justify the MHD approach occur if 
the turbulent fluctuations are small compared to the mean magnetic field, and having length- 
scales parallel to the mean field much larger than perpendicular length-scales. Treating wave 
frequencies that are low compared to the ion cyclotron frequency we enter the domain of 
"gyrokinetic approximation" which is commonly used in fusion plasmas e.g. [50,69], for which 
at length-scales larger than the ion gyroradius pi the incompressible shear-Alfven wave modes 
get decoupled from the compressive modes and can be described by the simple "reduced MHD" 
(RMHD) equations. These Alfven modes are most important for fast magnetic reconnection, what 
we discuss later in the review. 

In short, our considerations above confirm the generally accepted notion that the MHD 
approximation is adequate for most astrophysical turbulent plasmas at sufficiently large scales. 
In particular, the Goldreich-Srindhar [19] (henceforth GS95) theory of Alfvenic turbulence should 
be true for describing Alfvenic part of the MHD turbulent cascade^. For Alfvenic turbulence 
the eddies are elongated along magnetic field with the relation between the parallel and 
perpendicular dimensions due to the critical balance condition, namely, 

( 2 . 2 ) 

where Sug is the eddy velocity, while and £±_ are eddy scales parallel and perpendicular to 
the local direction of magnetic field, respectively. The notion of local magnetic field is the essential 
part of the modern understanding of Alfvenic turbulence and it was added to the GS95 picture by 
the later studies (LV99, [74,75]). The use of local magnetic field is expected as at small scale eddies 
can be influenced only by the magnetic field around them and not by the global mean field. 

A description of MHD turbulence that incorporates both weak and strong regimes was 
presented in LV99. In the range of length-scales where turbulence is strong, this theory implies 
that 

(2.3) 

Sui^UL(^^y (2.4) 

when the turbulence is driven isotropically on a scale Li with an amplitude uj^. As we see further, 
driving of turbulence by reconnection may be different from the isotropic driving assumed for the 
derivation of the expressions above. 

We do not discuss theories of Alfvenic turbulence that were develop to obtain the spectral 
index of —3/2 which was suggested by limited-resolution numerical simulations, e.g. in [75]^. 

The additional physics that was considered included, e.g. dynamical alignment [78], 
polarization intermittency [79], turbulence non-locality [80]. In particular [78] study predicts the 
the Kraichnan index of —3/2 [81,82] rather than Kolmogorov index —5/3 that follows from GS95. 
We feel that more recent high resolution numerical simulations (see [83,84]) provide results in 
agreement with the GS95 expectations, while the more shallow spectra are likely to be due to 
the bottleneck effect arising from MHD turbulence being less local compared to hydrodynamic 
one [85]. 

^We will concentrate on Alfvenic modes, while disregarding the slow and fast magnetosonic modes of MHD turbulence 
[70-72], which is possible as the backreaction of fast and slow inodes on Alfvenic cascade is insignificant [19,70,73]. 

^Low resolution numerical simulations are notorious in being ambiguous in terms of spectral slope. For instance, the initial 
compressible simulations suggested the spectral index of high Mach number hydrodynamic turbulence to be —5/3, which 
prompted theoretical attempts to explain this, e.g. [76]. However, further high resolution research [77] revealed that the 
flattering of the spectrum observed was the result of a bottleneck effect, which is more extended in compressible than in 
incompressible fluids. In the MHD simulations that are indicative of —3/2 spectrum, similar to the aforementioned low 
resolution hydrodynamic simulations showing —5/3 no bottleneck effect is seen. As the bottleneck is a physical effect, the 
fact that it is not seen in simulations to our mind means that it is just extended and higher resolution simulations are necessary. 
Therefore, choosing between theories on the basis of just spectral slope of low resolution simulations may be tricky. 
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The measurement in the solar wind show evidence for the —3/2 spectrum at the lAU from the 
Sun and —5/3 spectrum at distances larger than lAU [86]. We believe that the more relevant to 
MHD turbulence is the spectrum measured at larger distances where the there is less influence 
from the imbalance as well as the transient processes of spectrum evolution. While the discussion 
of the exact scaling of MHD turbulence is ongoing (see papers and comments by Perez et al. 
[87,88] and answers to them in [83,84]), we would like to stress that our results on recoimection 
marginally depend on the exact spectral index of turbulence. In LV99, which was developed when 
GS95 theory was far from being accepted, in the Appendix the reconnection rates were provided 
for arbitrary spectral indexes of turbulence and scale dependent anisotropies. 

More discussions of astrophysical turbulence can be found in recent reviews, e.g. [89-91]. 
In particular, there many additional effects are discussed, e.g. compressibility, effect of partial 
ionization as well as the effect of imbalance of turbulence. The latter may be a consequence of 
having sources and sinks of turbulent energy that are not coincident in space. All these effects are 
not of principal importance for our discussion of turbulent reconnection and therefore we do not 
discuss them here. 

Finally, we point out that we concentrate our attention on subAlfvenic turbulence as the 
reconnection of weakly perturbed magnetic fields is the natural generalization of the classical 
formulation of the reconnection problem. The opposite extreme is the turbulence in the 
dynamically unimportant magnetic field, where the magnetic field are reversing at the resistive 
dissipative scale. This is a degenerate example employed in the model of kinetic dynamo and it is 
of no interest for the reconnection research. If turbulence is superAlfvenic, magnetic field becomes 
dynamically important and stiff at the scale of LiM^^ [92] and the reconnection ideas below can 
be applied to such fields. 

The most important points of this section are 



• astrophysical fluids are generically turbulent, 

• MHD description is of Alfvenic turbulence is valid at sufficiently large scales, 

• we have an adequate theory of Alfvenic turbulence. 


In what follows we refer to these points dealing with the problem of turbulent reconnection. 


3. Turbulent Reconnection Model 

The model of turbulent reconnection in LV99 generalizes the classical Sweet-Parker model [93, 
94]'^. In the latter model two regions with uniform laminar magnetic fields are separated by thin 
current sheet. The speed of reconnection is given roughly by the resistivity divided by the sheet 
thickness, i.e. 

Vrecl-VlA. (3.1) 

For steady state reconnection the plasma in the current sheet must be ejected from the edge of the 
current sheet at the Alfven speed, Va- Thus the reconnection speed is 

Vrec2~VAAILoo, (3.2) 

where Lx is the length of the current sheet, which requires A to be large for a large reconnection 
speed. As a result, the overall reconnection speed is reduced from the Alfven speed by the square 
root of the Lundquist number, S = LxVA/Vr he. 

Vrec,SP = VAS-^/^. (3.3) 

The corresponding Sweet-Parker reconnection speed is negligible in astrophysical conditions as 
S may be 10^® or larger. 

The corresponding model of magnetic reconnection is illustrated by Figure 2. 

'^The basic idea of the model was first discussed by Sweet and the corresponding paper by Parker refers to the model as 
"Sweet model". 
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blow up 


Figure 2. Upper plot'. Sweet-Parker model of reconnection. The outflow is limited to a thin width <5, which is determined 
by Ohmic diffusivity. The other scale is an astrophysical scale 3> (5. Magnetic field lines are laminar. Middle plot: 
Turbulent reconnection model that accounts for the stochasticity of magnetic field lines. The stochasticity introduced by 
turbulence is weak and the direction of the mean field is clearly defined. The outflow is limited by macroscopic field line 
wandering. Low plot: An individual small scale reconnection region. [95]. 


Similar to the Sweet-Parker model, the LV99 model deals with a generic configuration, which 
should arise naturally as magnetic flux tubes try to make their way one through another. 
However, in the LV99 model the large-scale magnetic field wandering determines the thickness 
of outflow. Thus LV99 model does not depend on resistivity and can provide both fast and slow 
reconnection rates depending on the level of turbulence. 

To obtain the reconnection rate in LV99 model one should use the scaling relations for Alfvenic 
turbulence from § 2. A bundle of field lines confined within a region of width y at some particular 
point spreads out perpendicular to the mean magnetic field direcfion as one moves in either 
direction following the local magnetic field lines. The rate of field line diffusion is given by 


d{y^) . jy") ^ 

dx A|| 


(3.4) 


where Ajj"^ «is the parallel scale and the corresponding transversal scale, l±, is ~ 

and X is the distance along an axis parallel to the magnetic field. Therefore, using equation (2.3) 

one gets 


d{y^) 


U 




(3.5) 


where we have substituted for This expression for fhe diffusion coefficient will only 

apply when y is small enough for us to use the strong turbulence scaling relations, or in other 
words when (y^) < if (mj^/Va)^. 
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When the turbulence injection scale is less than the extent of the reconnection layer, i.e. Lx S> 
Li magnetic field wandering obeys the usual random walk scaling with Lx!Li steps and the 
mean squared displacement per step equal to L\{ulIVa)‘^. Therefore 

{yy^^~{L,x)^/\uLlVAf x>L, (3.6) 


Combining Eqs. (3.5) and (3.6) one can derive the thickness of the outflow A and obtain (LV99): 


Vrec ~ Va min 


1/2 


1 / 2 ' 


Ml, 


(3.7) 


where VaMI is proportional to the turbulent eddy speed. This reconnection rate represents a 
large fraction of the Alfven speed when Li and Lx are not too different and Ma is not too small. 

Due to the importance of the turbulent reconnection, it is advantageous to consider re¬ 
deriving the reconnection rates in another way. This was is based on the Lagrangian properties of 
magnetized plasma, in particular on the Richardson dispersion (see [96] and references therein). 

Richardson diffusion/dispersion can be illusfrafed with a simple hydrodynamic model. 
Consider the growth of the separation between two particles dl{t)/dt~ v(l), which for 
Kolmogorov turbulence is ~ where at is proportional to the energy cascading rate, i.e. 

at ~ V^/L for turbulence injected with superAlvenic velocity Vl at the scale L. The solution of 
this equation is 

= + at{t-to)fl (3.8) 


which at late times leads to Richardson diffusion/dispersion or ~ compared with l'^ 
for ordinary diffusion. Both terms "diffusion and dispersion" can be used interchangeably, but 
keeping in mind that the Richardson process results in superdiffusion (see [97] and references 
therein) we feel fhat it is advantageous to use the term "dispersion". Although the Richardson 
dispersion process was introduced for hydrodynamic turbulence a similar process is valid for 
magnetized fluids. We will not distinguish the magnetized and not magnetized case by name and 
instead of magnefic Richardson dispersion will use just Richardson dispersion. In magnetized 
turbulence Richardson dispersion is important in terms of spreading magnetic fields which 
provides a way to re-derive the LV99 relations. 

The fact that time dependence of the magnetic field diffusion induces magnefic reconnection, 
can be illustrated with the Sweet-Parker reconnection. There magnetic field lines are subjecf to 
Ohmic diffusion. The latter induces the mean-square distance across the reconnection layer that a 
magnetic field-line can diffuse by resistivity in a time t given by 

{y\t))~\t. (3.9) 

where A = <? jA-K a is the magnetic diffusivity. The field lines are advected out of the sides of the 
reconnection layer of length Lx at a velocity of order Va ■ Therefore, the time that the lines can 
spend in the resistive layer is the Alfven crossing time tA~ Lx/V a- Thus, field lines that can 
reconnect are separated by a distance 

^ {y^{tA)) ~ = Lx/Vs, (3.10) 

where S is Lundquist number. Combining Eqs. (3.2) and (3.10) one gets again the well-known 
Sweet-Parker result, Vrec = Va/^TS- 

The difference with the turbulent case is that instead of Ohmic diffusion one should use the 
Richardson one [32[. In this case the mean squared separation of particles (I®! (t) — 2:2 (t) 
where t is time, e is the energy cascading rate and (...) denote an ensemble averaging (see [98[). 
Por subAlfvenic turbulence e si u\/{VALi) (see LV99) and therefore analogously to Eq. (3.10) one 
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(3.11) 

where it is assumed that L<Li. Combining Eqs. (3.2) and (3.11) one obtains 

Vrec,LV99-VAiLIL,)^^^Ml. (3.12) 

in the limit of L < Li. Similar considerations allow to recover the LV99 expression for L > Li, 
which differs from Eq. (3.12) by the change of the power 1/2 to —1/2 and recover Eq. (3.7). 

4. Extending LV99 Reconnection Theory 

(a) Recent theoretical developments: rigorous mathematical approach 

Recently the LV99 notion of magnetic line wandering has played a central role in the extension 
of "general magnetic reconnection" (GMR) theory to turbulent plasmas. Recall that GMR theory 
[99,100] attempts to quantify the changes of magnetic connections between plasma elements. It 
is assumed in the standard approach to GMR that such changes occur only in narrow, sparsely 
distributed current layers or "diffusion regions" of small total volume. This assumption is invalid 
for turbulent plasmas. By tracking along field-lines wandering in space, [35] has developed an 
extended version of GMR theory valid for both laminar and turbulent plasmas. 

The study in [35] provides a rigorous mathematical treatment of the motion of magnetic field 
lines in turbulent plasmas. The slip source vector which is defined as the ratio of the curl of 
the non ideal electric field in the Generalized Ohm's Law and the magnetic field strength was 
introduced and it was demonstrated that this vector gives the rate of development of slip velocity 
per unit arc length of field line. It diverges at magnetic nulls, unifying GMR with magnetic 
null-point reconnection. In a turbulent inertial range the curl becomes extremely large while 
the parallel component is tiny, so that line slippage occurs even while ideal MHD is accurate. 
This means that ideal MHD is valid for a turbulent inertial-range only in a weak sense which 
does not imply magnetic line freezing (see also section 7). By rigorous estimates of the terms 
in the Generalized Ohm's Law for an electron-ion plasma the paper shows that all of the non¬ 
ideal terms (from collisional resistivity. Hall field, electron pressure anisotropy, and electron 
inertia) are irrelevant compared with the effects of turbulence and large-scale reconnection is 
thus governed solely by ideal dynamics. It is encouraging that in terms of magnetic reconnection 
the results of this study correspond to LV99 model and thus provide more rigorous theoretical 
foundations for turbulent reconnection. The results for the slippage velocity in [35] are identical to 
the expression of the reconnection velocities in LV99. Together with the earlier discussed results 
on Richardson dispersion in magnetic turbulence, these provide new outlook on the nature of 
magnetic reconnection in turbulent fluids. 

(b) Effect of energy dissipation in the reconnection layer 

In LV99 expressions were derived assuming that only a small fraction of the energy stored in 
the magnetic field is lost during large-scale reconnection and the magnetic energy is instead 
converted nearly losslessly to kinetic energy of the outflow. This can only be true, however, 
when the Alfvenic Mach number Ma is small enough. If Ma becomes large, then 

energy dissipation in the reconnection layer becomes non-negligible and there is a reduction of 
the outflow velocity (see ELVll). Note that even if A4 a is initially small, reconnection may drive 
stronger turbulence (see section 4c) and increase the fluctuation velocities wl in the reconnection 
layer. This scenario may be relevant to post-CME reconnection, for example, where there is 
empirical evidence that the energy required to heat the plasma in the recormection layer ("current 
sheet") to the observed high temperatures is from energy cascade due to turbulence generated by 
the reconnection itself [101]. In addition, Va within the reconnection layer will be smaller than 
the upstream values, because of annihilation of the anti-parallel components, which will further 
increase the Alfvenic Mach number. 
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The effect of turbulent dissipation can be estimated from steady-state energy balance in the 
reconnection layer: 

2 '^out^ ~ 2 '^A'^renLx — eLx^, (4-1) 

where kinetic energy carried away in the outflow is balanced against magnetic energy transported 
into the layer minus the energy dissipated by turbulence. Here we estimate the turbulent 
dissipation using the formula e = u^/VaLi for sub-Alfvenic turbulence [82]. Dividing (4.1) by 
A = LxVrec/vout, we get 

vLt = Vivout-2^^, (4.2) 

which is a cubic polynomial for Vout- The solutions are easiest to obtain by introducing the ratios 
/ = VoutIVA and r = 2M\{LxlLi) which measure, respectively, the outflow speed as a fraction 
of Va and the energy dissipated by turbulence in units of the available magnetic energy, giving 

r = f-f. (4.3) 

When r — 0, the only solution of (4.3) with / > 0 is / = 1, recovering the LV99 estimate Vout = 
Va for Ma <C 1. For somewhat larger values of r, / ~ 1 — {r 12), in agreement with the formula 
/ = (1 — r)^/^ that follows from Eq. (65) in ELVll, implying a slight decrease in Vout compared 
with Va- Note that formula (4.3) cannot be used to determine / for too large r, because it has 
then no positive, real solutions! This is easiest to see by considering the graph of r vs. /. The 
largest value of r for which a positive, real / exists is r-max = 2/v^ ~ 0.385 and then / takes 
on its minimum value fmin = 0.577. This implies that the LV99 approach is limited to 

M A sufficiently small, because of the energy dissipation inside the reconnection layer and the 
consequent reduction of the outflow velocity. This is not a very stringent limitation, however, 
because r is proportional to AtIf one assumes Lx ^ one may consider values oi A4 a up to 
0.662. Given the neglect of constants of order unity in the above estimate, we may say only that 
the LV99 approach is limited to At a F At the extreme limit of applicability of LV99, Vout is still 
a sizable fraction of Va, i.e. 0.577, not a drastically smaller value. 

The effect of the reduced outflow velocity may be, somewhat paradoxically, to increase the 
reconnection rate. The reason is that field-lines now spend a time Lxivout exiting from the 
reconnection layer, greater than assumed in LV99 by a factor of 1//. This implies a thicker 
reconnection layer A due to the longer time-interval of Richardson dispersion in the layer, greater 
than LV99 by a factor of (1//)^^^. The net reconnection speed Vrec = VoutA/Lx is thus larger 
by a factor of (1//)^^^. The increased width A more than offsets the reduced outflow velocity 
Vout- However, this effect can give only a very slight increase, at most by a factor of 3^~ 1.31 
for fmin = l/\/3. We see that for the entire regime Al^ < 1 where LV99 theory is applicable, 
energy dissipation in the reconnection layer implies only very modest corrections. It is worth 
emphasizing that "large-scale reconnection" in super-Alfvenic turbulence with At a > 1 is a very 
different phenomenon, because magnetic fields are then so weak that they are easily bent and 
twisted by the turbulence. Any large-scale flux tubes initially present will be diffused by the 
turbulence through a process much different than that considered by LV99. Eor a discussion of 
this regime, see [102]. 

(c) Reconnection in partially ionized gas 

On sufficiently small scales Alfvenic turbulence in the partially ionized gas is differs from our 
description provided in section 2. Due to viscosity caused by neutral atoms, the fluid viscosity 
becomes substantially larger than the fluid resistivity, which means that the Prandtl number of 
the fluid is high. Turbulence in high Prandtl number fluids has been studied numerically in [67, 
103,104] and theoretically in [95]. However, for our present discussion it is important that for 
scales larger than the viscous damping scale the turbulence follows the usual GS95 scaling and 
the considerations about Richardson dispersion and magnetic reconnection that accompany are 
valid at these scales. 
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In high Prandtl number media the GS95-type turbulent motions decay at the scale l_L^critr 
which is much larger than the scale of Ohmic dissipation. Thus over a range of scales less than 
l'±,crit to some much smaller scale magnetic field lines preserve their identity. To establish the 
range of scales at which magnetic fields perform Richardson diffusion one can observe that the 
transition to the Richardson dispersion is expected to happen when the field line separation 
reaches the perpendicular scale of the critically damped eddies l±^crit- The separation in the 
perpendicular direction starts with the scale rinit following the Lyapunov exponential growth 
with the distance I measured along the magnetic field lines, i.e. r^nit exp(i/Z|| crit)^ where 
corresponds to critically damped eddies with l±^crit- h seems natural to associate Vintt with the 
separation of the field lines arising from the action of Ohmic resistivity on the scale of the critically 
damped eddies 

^init ,crit /Ta 5 (4-4) 


where 77 is the Ohmic resistivity coefficient. 

At scales smaller than l±^crit fhe magnetic line separation obeys the laws established by 
Rechester & Rosenbluth [105]. The distance along the local magnetic field field over which 
anisotropic turbulence separates the magnetic field lines by l±^crit is the Rechester-Rosenbluth 
length (see [92]): 

^RR ^ ^W^crit (4-5) 

Taking into account Eq. (4.4) and that 

^^11 ,crzt /Ta 5 (4-6) 

where u is the viscosity coefficient, one can rewrite Eq. (4.5) as: 

Lrr~ Inherit Pi (4-7) 

where Pt = v/r] is the Prandtl number. This means that when the current sheets are much longer 
than Ljifj, then magnetic field lines undergo Richardson dispersion and according to [32] the 
recormection follows the laws established in LV99. At the same time on scales less than Lfiji 
magnetic reconnection maybe slow®. 


(d) Self-sustained Turbulent Magnetic Reconnection 

Recormection releases energy and induces outflows. Even if the initial magnetic field 
configuration is laminar, magnetic reconnection ought to induce turbulence due to the outflow 
(LV99, [106]). This effect was confirmed by observing the development of turbulence both in 
recent 3D Particle in Cell (PIC) simulations [107] and 3D MHD simulations [108,109]. 

In terms of MHD simulations, Beresnyak [108] was the first to study turbulent recormection 
with turbulence arising from the recormection itself. However, the periodic boundary conditions 
adopted in [108] limited the time span over which magnetic recormection can be studied and 
therefore the simulations focus on the process of establishing reconnection. 

Analytical description of the results in the framework of LV99 model was adopted by 
Beresnyak ( [108], private communication). Below we provide our theoretical account of the 
results in [108] using our imderstanding of LV99 turbulent recormection. We obtain expressions 
which are different from those by [108]. 

The logic of the derivation below is straightforward. As the opposite magnetic fluxes enters in 
contact, the width of the recormection layer A grows. The rate at which this happens is limited by 

^Incidentally, this can explain the formation of density fluctuations on scales of thousands of AU, that are observed in the 
ISM. 
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the mixing rate induced by the eddies at the scale A, i.e. 

Va 

A dt A 


(4.8) 


with a factor g which takes into account possible inefficiency in the diffusion process. As Va is a 
parf of the turbulent cascade, i.e. the mean value of « J ${ki)dki, where 


and Cfc is a Kolmogorov constant, which for ordinary MHD turbulence is calculated in [110], but 
in our special case may be different. If the energy dissipation rate e were time-independent, then 
the layer width would be implied by Eqs. (4.8) and (4.9) to grow according to Richardson's law 
A^ ~ et^. However, in the transient regime considered, energy dissipation rate is evolving. If the 
y-component of the magnetic field is reconnecting and the cascade is strong, then the mean value 
of the dissipation rate e is 

e^pvlyliAIVA), (4.10) 

where /3 is another coefficient measuring the efficiency of conversion of mean magnetic energy 
into turbulent fluctuations. This coefficient can be obtained from numerical simulations. 

The ability of the cascade to be strong from fhe very beginning follows from the large 
perturbations of the magnetic fields by magnetic reconnection, while magnetic energy can still 
dominate the kinetic energy. The latter factor that can be experimentally measured is given by a 
parameter r^. With this factor and making use of Eqs. (4.9) and (4.10), the expression for V/\ can 
be rewritfen in fhe following way: 

VA-CkrA^ylyVAfif'^ (4.11) 

where the dependences on fci ~ 1 jA cancel out. 

This provides the expression for fhe turbulent velocity at the injection scale Va 

Va « {CKTAf^VAyP^^^ (4.12) 

as a function of the experimentally measurable parameters of fhe system. Thus the growth of the 
turbulent reconnection zone is according to Eq.(4.8) 

^~gP^^\CKrAf^VAy (4.13) 

which predicts the nearly constant growth of the outflow region as seen in Fig.3 in [108]. 

For the steady state regime, one expects the outflow to play an important role. The equations 
for the reconnection rate were obtained in LV99 for the isotropic injection of energy. For the case 
of anisotropic energy injection of furbulence we should apply the following approach. Using Eq. 
(6.2) and identifying Va with the total velocity dispersion, which is similar to the use of Uots,turb 
in Eq. (6.1) one can get 

Vrec-VAiA/L:,)^/'^ (4.14) 


where the mass conservation condition provides the relation VrecLx ~ VAyA. Using the latter 
condition one gets 

Vrec~ VAyiCKTAf^^P (4.15) 

which somewhat slower than the rate at which the reconnection layer was growing initially. 


(e) Flares of Turbulent Reconnection 

On the basis of LV99 theory a simple quantitative model of flares was presented in [106]. There it is 
assumed that since stochastic reconnection is expected to proceed unevenly, with large variations 
in the thickness of the current sheet, one can expect that some unknown fraction of fhis energy 
will be deposifed inhomogeneously, generating waves and adding energy fo the local turbulent 
cascade. 
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For the sake of simplicity, the plasma density is assumed to be uniform so that the Alfven speed 
and the magnetic field sfrength are interchangeable. The nonlinear dissipation rate for waves is 


'^nonlinear ~ 


T liJ 2 , 

^ \ywave 7 2 T7 r 

«||kA 


(4.16) 


where the first rate is the self-interaction rate for the waves and the second is the dissipation rate 
induced by the ambient turbulence (see [111]). The important point here is that for the waves 
falls somewhere in the inertial range of fhe strong turbulence. Eddies at that wavenumber will 
disrupt the waves in one eddy turnover time, which is necessarily less than LjVA- Therefore, the 
bulk of the wave energy will go into the turbulent cascade before escaping from the reconnection 
zone. 

An additional simplification is achieved by assuming that some fraction e of the energy 
liberated by stochastic reconnection is fed into the local turbulent cascade. The evolution of the 
turbulent energy density per area is 

^ (^AV^')=tVlVrec-V^A^, (4.17) 

where the loss term covers both the local dissipation of turbulent energy, and its advection out 
of the reconnection zone. Since Vrec ~ Vturb and A ~ Lx{VIVa), if is possible to rewrite this by 
defining t = tVA/Lx so fhat 

■^Ma « eMA - Ma- (4.18) 

dr 

If e is a constant then 

(4.19) 


This implies that the time during which reconnection rate rises to (^^^Va is comparable to the 
ejection time from the reconnection region (~ Lx/Va)- 

Within this toy model e is not defined. Its value can be constrained through observations. 
Given that reconnection events in the solar corona seem to be episodic, with longer periods of 
quiescence, this is suggestive that e is very small, for example, depends strongly on the ratio of 
the thickness of the current sheet to Lx- In particular, if it scales as Ma to some power greater 
than two then initial conditions dominate the early time evolution. 

Another route by which magnetic reconnection might be self-sustaining via turbulence 
injection would be in the context of a series of topological knots in the magnetic field, each of 
which is undergoing reconnection. For simplicity, one can assume that as each knot undergoes 
reconnection it releases a characteristic energy into a volume which has the same linear dimension 
as the distance to the next knot. The density of the energy input into this volume is roughly 
LVaV! Lx, where here t is defined as fhe efficiency with which the magnetic energy is transformed 
into turbulent energy. Thus one gets 


VlV v'^ 
Lx T/fe 


(4.20) 


where Efc is the distance between knots and v' is the turbulent velocity created by the reconnection 
of fhe first knot. This process will proceed explosively ii v' >V or 

VlLkt>V^Lx. (4.21) 


The condition above is easy to fulfill. The bulk motions created by reconnection can generate 
turbulence as they interact with their surrounding, so e should be of order unity. Moreover the 
length of any current sheet should be at most comparable to the distance to the nearest distinct 
magnetic knot. The implication is that each magnetic reconnection event will set off its neighbors, 
boosting their reconnection rates from V^, sef by the environment, to t^^'^VAiLf^j Lx)^^"^ (as long 
as this is less than Va)- The process will take a time comparable to the crossing time Lx/Vl to 
begin, but once initiated will propagate through the medium with a speed comparable to speed of 
reconnection in the individual knots. The net effect can be a kind of modified sandpile model for 
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Figure 3. Visualization of reconnection simulations in [33,34], Left ipaneh Magnetic field in the reconnection region. 
Central panel: Current intensity and magnetic field configuration during stochastic reconnection. The guide field is 
perpendicular to the page. The intensity and direction of the magnetic field is represented by the length and direction 
of the arrows. The color bar gives the intensity of the current. Right panel: Representation of the magnetic field in the 
reconnection zone with textures. 


magnetic reconnection in the solar corona and chromosphere. As the density of knots increases, 
and the energy available through magnetic reconnection increases, the chance of a successfully 
propagating reconnection front will increase. 

(f) Relativistic reconnection 

Magnetic turbulence in a number of asfrophysical highly magnetized objects, accretion disks 
near black holes, pulsars, gamma ray bursts may be in the relativistic regime when the Alfven 
velocity approaches that of light. The equations that govern magnetized fluid in this case look 
very different from the ordinary MHD equations. However, studies by [112] and [113] show 
that for bofh balanced and imbalanced turbulence, the turbulence spectrum and turbulence 
anisotropies are quite similar in this regime and the non-relativistic one. This suggests that the 
Richardson dispersion and related processes of LV99-fype magnetic reconnection should cary 
on to the relativistic case (see [114]). This prediction was confirmed by the recent numerical 
simulations Makoto Takomoto (2014, private communication) who with his relativistic code 
adopted the approach in [33] and showed that the rate of 3D relafivisfic magnetic reconnection 
gets independent of resistivity. 

The suggestion that LV99 is applicable to relativistic reconnection motivated the use of the 
model for explaining gamma ray bursfs in [115] and [36] studies and in accretion disks around 
black holes and pulsars studies [116,117]. Now, as the extension of the model to relativistic 
case has be confirmed these and other cases where the relativistic analog of LV99 process was 
discussed to be applicable (see [21]) are given numerical support. 

Naturally, more detailed studies of both relativistic MHD turbulence and relativistic magnetic 
reconnection are required. It is evident that in magnetically-dominated, low-viscous plasmas 
turbulence is a generic ingredient and thus it must be taken into account for relativistic magnetic 
reconnection. As we discuss elsewhere in the review the driving of turbulence may by external 
forcing or it can be driven by reconnection itself. 

5. Numerical Testing of Turbulent Reconnection Theory 

Figure 3 illustrates results of numerical simulations of turbulent reconnection with turbulence 
driven using wavelets in [33] and in real space in [34]. 

As we show below, simulations in [33,34] provided very good correspondence to the 
LV99 analytical predictions for the dependence on resistivity, i.e. no dependence on resistivity 
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Figure 4. Left Panel The dependence of the reconnection velocity on the injection power for different simuiations with 
different drivings. The predicted LV99 dependence is also shown. Pinj and ki„j are the injection power and scale, 
respectively, Bz is the guide field strength, and rju the value of unifor resistivity coefficient. Right Panel The dependence 
of the reconnection velocity on the injection scale. From [34], 
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1 /2 

for sufficiently strong turbulence driving, and the injection power, i.e. Vrec ~ ■ The 

corresponding dependence is shown in Figure 4, left panel. 

The simulations did not reveal any dependence on the strength of the guide field Bz (see 
Figure 4). To address this dependence, in the limit where the parallel wavelength of the strong 
turbulent eddies is less than the length of the current sheet, we can rewrite the reconnection speed 
as 


Vr, 


1 

VAx J fc||K4 


(5.1) 


Here F is the power in the strong turbulent cascade. Lx and Vax are the length scale and Alfven 
velocity in the direction of the reconnecting field, and Va is the total Alfven velocity, including 
the guide field. 

In a physically realistic situation, the dynamics that drive the turbulence, whatever they are, 
provide a characteristic frequency and input power. Since the guide field enters only in the 
combination feyVA, i.e. through the eddy turn over rate, this implies that varying the guide 
field will not change the reconnection speed. In the simulations the periodicity of the box in the 
direction of the guide field complicates the analysis (see more discussion in [118]). 

The injection of energy in LV99 is assumed to happen at a given scale and the inverse cascade 
is not considered in the theory. Therefore it is not unexpected that the measured dependence on 
the turbulence scale differs from the predictions. In fact, it is a bit more shallow compared to the 
LV99 predictions (see Figure 4, right panel). 

The left panel of Figure 5 shows the dependence of the reconnection rate on explicit uniform 
viscosity obtained from the isothermal simulations of the magnetic reconnection in the presence 
of turbulence [34]. The open symbols show the reconnection rate for the laminar case when there 
was no turbulence driving, while closed symbols correspond to the mean values of reconnection 
rate in the presence of saturated turbulence. All parameters in those models were kept the same, 
except the uniform viscosity, which varied from 10“^ to 10“^ in the code imits. We notice the lack 
of any scaling for the laminar case, which is somewhat in contradiction to the scaling Vrec oc 
derived in [119]. We should remind, that the authors introduced the viscosity dependency using 
the energy equation balance, which cannot be applied in the isothermal case. They also stress 
that the proper scaling might be sensitive to the chosen boundaries, which in theirs numerical 
tests where closed. In the models presented in Figure 5 we use outflow boundaries. The viscosity 
scaling for the case when turbulence is present is shown by closed symbols. This scaling is also 
Vrec OC , but Can be explained rather as the effect of the finite inertial range of turbulence 
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Figure 5. Left panel. The dependence of the reconnection velocity on uniform viscosity in the 3D isothermai models of 
Sweet-Parker reconnection (open symbols) and reconnection enhanced by the presence of turbulence (closed symbols) 
from [34], The dependence on viscosity is negligible in the laminar case, while in the presence of turbulence the 
reconnection rate exhibits the scaling See text for explanation. 

Right panel. The reconnection rate in models with anomalous resistivity for Sweet-Parker case 
(filled circles) and in fhe presence of furbulence (filled diamonds). We observe no dependence of 
the reconnection rate on the strength of anomalous effects. In the Sweet-Parker case the 
anomalous resistivity is not turned on, since the maximum current density is below the 
threashold of the anomalous resistivity model we used. See [33] for more detailed description. 


than the effect of energy balance affected by viscosity or boimdary conditions. For an extended 
range of motions, LV99 does not predict any viscosity dependence, if fhe dissipation scale lays 
much below the scale of current sheet. However, for numerical simulations the range of turbulent 
motions is very limited and any additional viscosity decreases the resulting velocity dispersion 
and therefore the field wandering fhus affecting fhe reconnection rafe. 

LV99 predicted that in the presence of sufficientiy sfrong turbulence, plasma effects should not 
play a role. The accepted way to simulate plasma effects within MHD code is to use anomalous 
resistivity. The results of the corresponding simulations are shown in the right panel of Figure 5 
and they confirm that the change of the anomalous resistivity does not change the reconnection 
rate. 

Within the derivation adopted in LV99 current sheet is broad with individual currents 
distributed widely within a three dimensional volum and the turbulence within the reconnection 
region is similar to the turbulence within a statistically homogeneous volume. Numerically, the 
structure of the reconnection region was analyzed by Vishniac et al. [120] based on the numerical 
work by Kowal et al. [33]. The results support LV99 assumptions with reconnection region being 
broad, the magnetic shear is more or less coincident with the outflow zone, and the turbulence 
within it is broadly similar to turbulence in a homogeneous system. In particular, this analysis 
showed that peaks in the current were distributed throughout the reconnection zone, and that the 
width of these peaks were not a strong function of fheir sfrength. The illustration of the results 
is shown in Figure 6 which shows histograms of magnetic field gradienfs in the simulations with 
strong and moderate driving power, with no magnetic field reversal but with driven turbulence, 
and with no driven turbulence at all, but a passive magnetic field reversal (i.e. Sweet-Parker 
reconnection). A few features stand out in this figure. Firsf, all the simulations with driven 
turbulence have a roughly gaussian distribution of magnetic field gradients. In the case with no 
field reversal (panel c) the peak is narrow and symmetric around zero. In the presence of a large 
scale field reversal the peak is slightly broadened, and skewed. It is turbulent reconnection does 
not produce any strong feature corresponding to a preferred value of fhe magnetic field gradienf. 
Instead one sees a systematic bias towards large positive values. We conclude from the lack of 
coherent features within the outflow zone, and the broad distribution of values of the gradient 
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Figure 6. These figures show histograms of the gradient of the reversing component of the large scale magnetic field in 
the direction normal to the unperturbed current sheet, i.e. dyB^. Upper left panel is for the highest power simulation, 
P=1. Upper right panel is for P=0.5. Lower left is for P=t but with no large scale magnetic field reversal, i.e. simply locally 
driven strong turbulence. Bins with twice the number of cells as the corresponding bin with the opposite sign of dyB^ are 
shown in green. Lower right shows the first simulation in the absence of turbulent forcing. From [120]. 


of the magnetic field, that the current sheet and the outflow zone are roughly coincident and this 
volume is filled with turbulent structures. 

As we discussed, the LV99 model is intrinsically related to the concept of Richardson 
dispersion in magnetized fluids. Thus by testing the Richardson diffusion of magnetic field, one 
also provides tests for the theory of turbulent reconnection. 

The first numerical tests of Richardson dispersion were related to magnetic field wandering 
predicted in LV99 [95,108,121]. In Figure 7 we show the results obtained in [95]. There we 
clearly see different regimes of magnetic field diffusion, including the y ~ regime. This is 
a manifestation of the spatial Richardson dispersion. 

As we discussed in section 3, the LV99 expressions can be obtained by applying the concept 
of Richardson dispersion to a magnetized layer. Thus by testing the Richardson diffusion of 
magnetic field, one also provides tests for the theory of turbulent reconnection. 

The numerical tests of Richardson dispersion in space correspond to magnetic field wandering 
predicted in LV99. In Figure 7 we show the results obtained in [95]. There we clearly see the 
Richardson regime corresponding to y ~ regime (see more discussion in ELVll). 

A direct testing of the temporal Richardson dispersion of magnetic field-lines was performed 
recently in [122]. For this experiment, stochastic fluid trajectories had to be tracked backward in 
time from a fixed point in order to determine which field lines at earlier times would arrive to that 
point and be resistively "glued together". Hence, many time frames of an MHD simulation were 
stored so that equations for the trajectories could be integrated backward. The results of this study 
are illustrated in Figure 7. The left panel shows the trajectories of the arriving magnetic field¬ 
lines, which are clearly widely dispersed backward in time, more resembling a spreading plume 
of smoke than a single "frozen-in" line. Quantitative results are presented in the right panel. 
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Figure 7. Left panel. Stochastic trajectories that arrive at a fixed point in the archived MHD fiow, coior-coded red, green, 
and biue from earlier to later times. From [122], Right panel. Mean-square dispersion of fieid-iines backwards in time, 
with red for direction parailei and blue for direction perpendicular to the local magnetic field. From [122]. 


which plots the root-mean-square line dispersion in directions both parallel and perpendicular to 
the local mean magnetic field. Times are in units of the resistive time 1/jrms determined by the 
rms current value and distances in units of the resistive length X/jrms- The dashed line shows 
the standard diffusive estimate 4At, while the solid line shows the Richardson-type diffusion, 
the power-law is a bit altered by the numerical effects® We would like to stress that whatever 
plasma mechanism of line-slippage holds at scales below the ion gyroradius— electron inertia, 
pressure anisotropy, etc.—will be accelerated and effectively replaced by the ideal MHD effect of 
Richardson dispersion. 

As we discussed in Section 4c the self-sustained turbulent reconnection where the turbulence 
is generated by the reconnection itself can be quantified using the predictions of the LV99 theory. 
Below we compare the prediction given by Eq. (4.15) against the results of recent simulations 
illustrated by Figure 8. The figure shows a few slices of the magnetic field strength \B\ through 
the three-dimensional computational domain with dimensions Lx = 1.0 and Ly = Lz= 0.25. The 
simulation was done with the resolution 2048 x 512 x 512. Open boundary conditions along the 
X and Y directions allowed studies of sfeady state turbulence. At the presented time t = 1.0 
the turbulence strength increased by two orders of magnifude from its initial value of ~ 
10~‘^Emag- Initially, only the seed velocity field at the smallest scales was imposed (a random 
velocity vector was set for each cell). We expect that most of the injected energy comes from the 
Kelvin-Helmholtz instability induced by the local interactions between the reconnection events, 
which dominates in the Z-direction, along which a weak guide field is imposed {Bz = O.lBx). As 
seen in the planes perpendicular to Bx in Figure 8, Kelvin-Helmholtz-like structures are already 
well developed at time t = 1.0. Turbulent structures are also observed within the XY-plane, which 
probably are generated by the strong interactions of the ejected plasma from fhe neighboring 
reconnection events. More detailed analysis of the spectra of turbulence and efficiency of the 
Kelvin-Helmholtz instability as the turbulent injection mechanism are presented in [109]. 

The Kelvin-Helmholtz instability due to the interactions of the outflows from neighboring 
reconnection events, which takes place in our simulations, is somewhat different from that in the 
current sheet of Sweet-Parker reconnection, which has been theoretically predicted in [123]. In 
the laminar reconnection, the profile of the outflow velocity has its maximum in the middle of the 


^Due to the bottleneck effect the measured magnetic energy spectrum is k [85] and this spectrum corresponds to 
Richardson dispersion dependence. 
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current sheet and quickly decays along the direction parallel to the reconnecting magnetic field 
component. This configuration creates naturally two shear layers in which the Kelvin-Helmholtz 
instability may develop if fhe oufflow velocity exceeds the Alfven speed associated with the 
upstream magnetic field. In order fo confirm fhe predictions obtained in [123] we would need 
simulations or observations of fhe thin current sheets with very large resolutions. 



Figure 8. Visualization of the model of turbulence generated by the seed reconnection from [f 09]. Three different outs 
(one XY plane at Z=-0.f and two YZ-planes at X=-0.25 and X=0.42) through the computational domain show the 
strength of magnetic field |.B| at the evolution time t = 1.0. Kelvin-Helmholtz-type structures are well seen in the planes 
perpendicular to the reconnecting magnetic component In the Z direction, the Kelvin-Helmholtz instability is slightly 
suppressed by the guide field of the strength = O.IB^ (with B^ = 1.0 initially). The initial current sheet is located 
along the XZ plane at Y=0.0. A weak ~ lO~'^Emag) random velocity field was imposed initialy in order to seed 
the reconnection. 


In Figure 9 we show the growfh of the turbulent region thickness for models ran with the 
same parameters, one in periodic domain, and another with open boundary conditions. In the 
case of periodic box, the reconnected flux and generated turbulence are accumulated in the region 
near current sheet. Since there is no outflow, the thickness of fhis region grows linearly with the 
estimated growth rate of about 0.026 (see blue line in Fig. 9). Once we allow the reconnected flux 
and furbulence to be ejected along the reconnecting magnetic field, the thickness of furbulenf 
region saturates after about 2.0 Alfven fime units at level of 0.025 L, where L is the longitudinal 
size of the domain. These values are in agreement with the estimates from Eqs. (4.13) and (4.15). 


6. Observational Testing of Turbulent Reconnection 
(a) Solar Reconnection 

To quantify solar reconnection one should accept that the energy is injected by reconnection 
and turbulence is driven by magnetic reconnection. In this situation one can expect substantial 
changes of the magnetic field direction corresponding to strong turbulence. Thus it is natural to 
identify fhe velocities measured during the reconnection events with the strong MHD turbulence 
regime. In other words, one can use: 


Vr, 


^obs ,turb{^inj / 


( 6 . 1 ) 
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Time 


Figure 9. The growth of the turbulent region width for two models similar to one presented in Fig. 8 but with periodic and 
open domain. From [109]. 


where Uobs^turb is the spectroscopically measured turbulent velocity dispersion. Similarly, the 
thickness of the reconnection layer should be defined as 

^ « Lx{Uobs,turb|VA){L^nJ|L^)^^^. ( 6 . 2 ) 

The expressions given by Eqs. (6.1) and (6.2) can be compared with observations in [124]. 
There, the widths of the reconnection regions were reported in the range from O.OSLx up to 0.16La; 
while the the observed Doppler velocities in the units of Va were of the order of 0.1. It is easy to 
see that these values are in a good agreement with the predictions given by Eq. (6.2)^ 

If we talk about unique predictions that radically differ from LV99 and the present day plasma 
reconnection models then the LV99 prediction of fhe triggering of reconnection by wave packets 
coming from the adjacent reconnectionsites should be singled out. Thus a particular series of 
solar observations is important. In [125], authors explaining quasi-periodic pulsations in observed 
flaring energy releases at an active region above the sunspot, proposed that the wave packets 
arising from the sunspots can trigger such pulsations. This is exactly what is expected within the 
LV99 model. 

The criterion for the application of LV99 theory is that the outflow region is much larger than 
the ion Larmor radius pi. This is definitely satisfied for the solar atmosphere where the ratio 
of A to pi can be larger than 10®. Plasma effects can play a role for small scale reconnection events 
within the layer, since the dissipation length based on Spitzer resistivity is ~ 1 cm, whereas pi ~ 
10® cm. However, as we discussed earlier, this does not change the overall dynamics of turbulent 
reconnection. 

(b) Solar Wind 

Reconnection throughout most of the heliosphere appears similar to that in the Sun. For example, 
there are now extensive observations of reconnection jets (outflows, exhausts) and strong current 
sheets in the solar wind [126]. The most intense current sheets observed in the solar wind are very 
often not observed to be associated with strong (Alfvenic) outflows and have widths at most a 
few tens of the proton inertial length 5i or proton gyroradius pi (whichever is larger). Small-scale 
current sheets of this sort that do exhibit observable reconnection have exhausts with widths 

we associate the observed velocities with isotropic driving of turbulence, which is unrealistic for the present situation, then 
a discrepancy with Eq. (6.2) would appear. Because of that [124] did not get quite as good quantitative agreement between 
observations and theory as we did, but still within observational uncertainties. 
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at most a few hundreds of ion inerfial lengths and frequently have small shear angles (strong 
guide fields) [127,128]. Such small-scale reconnection in the solar wind requires collisionless 
physics for its description, but the observations are exactly what would be expected of small- 
scale reconnection in MHD turbulence of a collisionless plasma [129]. Indeed, LV99 predicted that 
the small-scale reconnection in MHD turbulence should be similar to large-scale reconnection, 
but with nearly parallel magnetic field lines and with "outflows" of fhe same order as the local, 
shear-Alfvenic turbulent eddy motions. It is worth emphasizing that reconnection in the sense 
of flux-freezing violation and disconnection of plasma and magnetic fields is required at every 
point in a turbulent flow, not only near the most intense current sheets. Otherwise fluid motions 
would be halted by the turbulent tangling of frozen-in magnetic field lines. However, except 
at sporadic strong current sheets, this ubiquitous small-scale turbulent reconnection has none 
of the observable characteristics usually attributed to reconnection, e.g. exhausts stronger than 
background velocities, and would be overlooked in observational studies which focus on such 
features alone. 

However, there is also a prevalence of very large-scale reconnection events in the solar wind, 
often associated with interplanetary coronal mass ejections and magnetic clouds or occasionally 
magnetic disconnection events at the heliospheric current sheet [126,130]. These events have 
reconnection outflows with widths up to nearly 10® of the ion inertial length and appear to be in 
a prolonged, quasi-stationary regime with reconnection lasting for several hours. Such large-scale 
reconnection is as predicted by the LV99 theory when very large flux-structures with oppositely- 
directed components of magnetic field impinge upon each other in the turbulent environment of 
the solar wind. The "current sheet" producing such large-scale reconnection in the LV99 theory 
contains itself many ion-scale, infense current sheets embedded in a diffuse furbulenf background 
of weaker (but still substantial) current. Observational efforts addressed to proving/disproving 
the LV99 theory should note that it is this broad zone of more diffuse currenf, not the sporadic 
strong sheets, which is responsible for large-scale furbulent reconnection. Note that the study 
[122] showed that standard magnetic flux-freezing is violated at general points in turbulent MHD, 
not just at the most intense, sparsely distributed sheets. Thus, large-scale reconnection in the solar 
wind is a very promising area for LV99. 

Preliminary comparisons between such events in MHD turbulence and in the high-speed solar 
wind have yielded very promising results [131]. Criteria can be employed that are designed 
specifically to look for large-scale reconnection. For example, the "partial-variance of increments" 
(PVI) criterion recently proposed by [132] can be adapted for this purpose, by considering 
magnetic increments over inertial-range separation distances rather than ion-scale distances and, 
possibly also, with coarse-graining of fhe magnetic field to eliminate smaller-scale features. A 
similar modification may be made to the criterion of Gosling [126], which identifies reconnection 
events by roughly Alfenic-jettrng plasma bounded on one side by correlated changes in the 
antiparallel components of u and B and by anti-correlated changes in those components on the 
other side. Here the criterion may be modified by requiring that the two large changes must be 
separated spatially by inertial-range lengths, i.e. essentially by conditioning on a broad outflow 
jet. 

Examples of some events yielded by this latter Gosling-type criterion are shown in Figure 10. 
The top panels of the figures show a typical event selected from the JHU turbulence database, 
which archives the output of a 1024® pseudo-spectral simulation of the incompressible MHD 
equations. The bottom panel shows a similar event obtained from a study of a fast solar wind 
stream, 2008 January 14 04:40:00 - January 21 03:20:00, using three-second resolution Wind 
spacecraft observations from the Magnetic Field Investigation (MFI) and 3D Plasma Analyzer 
(3DP) experiments. The left panels show magnetic field componenfs and the right panels show 
velocity components, both rotated into the local mrnimum-variance-frame [133] plotted versus 
space for ID cuts through the MHD simulation and versus time for the spacecraft data. The JHU 
MHD data are in the arbitrary units of the simulation, for which the rms magnetic field strength 
b' = 0.24, the magnetic integral length Lf, = 0.35, and the resistive dissipation length rji, — 0.0028. 
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The units for the Wind data are nanotesla (nT) for the magnetic field, kilometer per second (km/s) 
for velocify, and minutes for time. Average solar wind conditions were speed u= 660 km/s, 
magnetic field sfrength B = 4.4 nT, proton number density rip = 2.4 cm“^, Alfven speed Va =62 
km/s, and proton beta /3p = 1.2. The outer scale of the turbulent inertial-range (boundary with 
the 1 // spectral range) is 33 mins and the inner scale (a few ion gyroradii) around 10 s. 

The event from the MHD database was found by searching for "Gosling events" that show 
opposing changes in u and B within a distance of 0.196, about half an integral length. The event 
from the high-speed solar wind was found by applying the same criterion for separation of 400 s. 
Interestingly, neither of these events show the "double-step" structure, with an intermediate 
plateau of reversing magnetic field component, which often characterize the events identified by 
Gosling [126], although other events we have foimd do show this structure. Most importantly, 
both events show the features expected of large-scale reconnection, with a sizable magnetic 
reversal over an inertial-range length and with a corresponding outflow in the same direction 
and of the same width. This makes both events likely candidates for turbulent reconnection. In 
the case of the MHD database event, this interpretation can be verified from the simulation data. 
A detailed study in preparation (Eyirrk et al., in prep.) shows that the MHD event presented 
in the top panels of Figure 10 accords well with the predictions of the LV99 theory and has 
the expected morphological features: a wide (inertial-range scale) outflow jet, a distribution of 
small-scale current sheets rather than a single dominant sheet, turbulent wandering of magnetic 
field-lines, and Richardson dispersion of field-lines normal to the reversal direction. It is therefore 
natural to identify the similar events in the solar wind as turbulent reconnection as well. This 
identification is strengthened by the similar statistical rates of occurrence of such events at 
corresponding scales, as observed also in previous studies of inertial-range magnetic increments 
in MHD turbulence and the solar wind [134]. The high-speed solar wind is presumably full of 
such turbulent reconnection events, across its broad spectrum of inertial-range length-scales. 

We note, that the situation for applicability of LV99 generally gets better with increasing 
distance from the sun, because of fhe great increase in scales. For example, reconnecting flux 
structures in the inner heliosheath could have sizes up to ~100 AU, much larger than the ion 
cyclotron radius ~ 10^ km [135]. 

(c) Parker spiral and Heliospheric Current Sheet 

More recently, ]35] discussed some implications of LV99 for heliospheric reconnection, in 
particular for deviations from the Parker spiral model of interplanefary magnetic field. Note, 
that the [136] spiral model of the interplanetary magnetic field, which is one of the most famous 
applications in asfrophysics and space science of the "frozen-in" principle for magnetic field 
lines. The model has been shown to be approximately valid when taking into account solar cycle 
variations in source magnetic field strength and latitude/time variation in solar wind speeds. 
Nevertheless, [136] concluded his paper with a "warning to the reader against taking too literally 
any of the smooth idealized models which we have constructed in this paper". 

]137] had studied the magnetic geometry and found "notable deviations" from the spiral 
model. [137] studied daily averages of magnetic field observations of Voyager 1 and 2 in the 
ecliptic plane at solar distances R=l-5 AU during a period of increasing solar activity in the 
years 1977-1979. In contrast to the Parker predictions for radial magnetic field component radial 
dependencies Bfi ~ R~^ and azimuthal component B^ ~ R~^, [137] found Bji ~ and 

Bt ~ ■ [137] attributed the observed deviations to "temporal variations associated with 

increasing solar activity, and to the effects of fluctuations of fhe field in the radial direction". These 
early observations were recently confirmed by [138], who presented evidence on the breakdown 
of the Parker spiral model for time- and space-averaged values of fhe magnetic field from several 
spacecraft (Helios 2, Pioneer Venus Orbiter, IMPS, Voyager 1) in the inner heliosphere at solar 
distances 0.3-5 AU and in the years 1976-1979. ]138] interpret their observations as due to "a 
quasi-continuous magnetic reconnection, occurring both at the heliospheric current sheet and at 
local current sheets inside the IMF sectors". They present extensive evidence that most nulls of 
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Figure 10. Candidate Events for Turbulent Reconnection. 

MHD Turbulence Simulation (Top Panels) and High-Speed Solar Wind (Bottom Panels). The left panels 
show magnetic field components and the right panels velocity components, both rotated into a 
local minimum-variance frame of the magnetic field. The component of maximum variance in 
red is the apparent reconnecting component, the component of medium variance in green is the 
nominal guide-field direction, and the minimum-variance direction in blue is perpendicular to 

the reconnection layer. 


BR and BT , where reconnection may occur, are not associated to the heliospheric current sheet. 
They as well observe a rapid disappearance of the regular sector structure at distances past 1 
AU, which they attribute to £)turbulent processes in the inner heliosphere. [35] estimated the 
magnetic field slippage velocities and related the deviation from Parker original predictions to 
LV99 reconnection. 

In addition, [35] analyzed the data relevant to the region associated with the broadened 
heliospheric current sheet (HCS), noticed its turbulent nature and provided arguments on 
the applicability of LV99 magnetic reconnection model to HCS. This seems to be a very 
promising direction of research to study turbulent reconnection in action using in situ spacecraft 
measurements. 


7. Implications 
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(a) Magnetic flux freezing in the presence of turbulent reconnection 

The concept of flux freezing was first proposed by Hannes Alfven in 1942, and the principle 
of frozen-in field lines has provided a powerful heuristic [2,139]. However, if strictly valid, flux 
would forbid magnetic reconnection, because field-lines frozen into a continuous plasma flow 
cannot change their topology. If magnefic reconnection were a phenomenon isolated to to regions 
of special magnetic flux topology or other special conditions, then it would be possible to use 
flux freezing for generic magnetic field conditions. However, LV99 model suggests that magnetic 
reconnection happens everywhere in magnetized turbulent fluids. This means the ubiquitous 
violation of flux freezing in magnetized turbulence. 

Standard mathematical proofs of flux-freezing in MHD always assume, implicitly, that velocity 
and magnetic fields remain smoofh as ?7 —>■ 0. However, MHD solutions with small resistivities 
and viscosities (high magnetic and kinetic Reynolds numbers) are generally turbulent. These 
solutions exhibit long ranges of power-law spectra corresponding to very non-smooth or "rough" 
magnetic and velocity fields. Fluid particle (Lagrangian) trajectories in such rough flows are 
known to be non-unique and stochastic (see [140-145], and, for reviews, [98] and [146]). It view 
of fhe above, it is immediately clear as a consequence that standard flux-freezing cannot hold in 
turbulent plasma flows. After all, the usual idea is that magnetic field-lines at high conductivity 
are tied to the plasma flow and follow the fluid motion. However, if the latter is non-imique and 
stochastic, then which fluid element will the field-line follow? 

For a laminar velocity field, this diffusion effect is small. It is not hard to see that a pair of field 
lines will attain a displacement r{t) apart imder the combined effect of advection and diffusion 
obeying 

^(r^) = 12A-F2(r-<5u(r)) 

where 5u(r) is the relative advection velocity at separation r. Thus, 

^(r^)<12A-F2||Vu||(r^), 

where || Vu|| is the maximum value of the velocity-gradient Vu. It follows that two lines initially 
at the same point, by time t can have separated at most 


(r^(t)) <6A 


g2||Vu||t _ ^ 

liVull 


(7.1) 


If we thus consider a smooth laminar flow with a fixed, finite value of || Vu||, then {r'^(t)) —>■ 0 as 
A —>• 0. Under such an assumption, magnetic field lines do not diffuse a far disfance away from 
the solution of the deterministic ordinary differential equation dx/dt — u(x, t), and the magnetic 
line-diffusion becomes a negligible effect. In that case, magnetic flux is conserved better as A 
decreases. 

However, in a turbulent flow, the above argument fails! The inequality (7.1) still holds, of 
course, but it no longer restricts the dispersion of field-lines under the joint action of resistivity 
and advection. As is well-known, a longer and longer inertial range of power-law spectrum 
E{k) occurs as viscosity i/ decreases and the maximum velocity gradient ||Vu|| becomes larger 
and larger. In fact, energy dissipation e = i^|j Vu||^ is observed to be non-vanishing as t 0 in 
turbulent flow, requiring velocity gradients to grow unboimdedly. Estimating || Vu|| ~ , 

the upper bound (7.1) becomes 


< 6A(z.//e)^^^[exp(2t(£/i/)^^^) — 1]. 


(7.2) 


This bound allows unlimited diffusion of field-lines. Consider first the case X = u or Pt = 1, for 
simplicity, where Richardson's theory implies that 

{r‘^{t))r^l2Xt + et^. (7.3) 


The rigorous upper bound always lies strictly above Richardson's prediction and, in fact, goes 
to infinity as i/ = A —>■ 0! The case of large Prandtl number is jusf slightly more complicated, as 
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previously discussed in §(c). When Pt^ 1, the inequality (7.2) holds as an equality for times 

t < ttrans with 


In(Pt) 

ttrans 2{e / u)' 

This is then followed by a Richardson despersion regime 

(r^(t)) ~ + e{t - Uransf, t > ttrans 


(7.4) 

(7.5) 


assuming that the kinetic Reynolds number is also large and a Kolmogorov inertial range exists at 
scales greater than the Kolmogorov length Once again, the upper bound (7.2) is much 

larger than Richardson's prediction and, at times longer than ttrans, the dispersion of field lines 
is independent of resistivity. 


(b) Making Goldreich-Sridhar model self-consistent 

Historically a lot of reconnection research was aimed to obtain the Holy Grail number of 
reconnection speed, which on the basis of solar flare observations was determined to be O.IVa. 
This recormection speed has been recently claimed to be attained in a number of plasma 
simulations (see [16]). We claim, however, that to make any model of strong turbulence self- 
consistent the velocity of O.IVa is insufficient. Below we show this for the GS95 model by 
reproducing the arguments in LV99. Magnetic reconnection is required for free mixing of 
magnetic field lines, which is a part of the GS95 picture of turbulence. In fact, the critical balance 
that is the corner stone of the GS95 model can be derived from the equality of times for mixing 
of magnetic field lines perpendicular to the local direction of magnetic field and the period of 
the Alfven wave that this mixing induces. Therefore we consider magnetic reconnection within 
magnetic eddies elongated along the local magnetic field direction. 

It is possible to see that within the GS95 picture the reconnection happens with nearly parallel 
lines with magnetic pressure gradient being reduced by a factor , since only reversing 

component is available for driving the outflow. At the same time the length of the contracted 
magnetic field lines is also reduced from by Therefore the acceleration is '’'ejtct^'±/h\' 

As a result, the Newtons' law gives yA^±/^\\^'^eject^±/h\' provides the result for the 
ejection rate Va/1\\ - The length over which the magnetic eddies intersect is l± and the 

rate of reconnection is Vrecll±- For the stationary reconnection this gives Vrec ~ VaI-±/1\\/ which 
provides the reconnection rate Va/1\\ - The latter rate is exactly the rate of the eddy turnovers in 
GS95 turbulence, which shows that it is fast magnetic reconnection that makes the GS95 picture 
self-consistent. In the case of trans-Alfvenic turbulence this means that the reconnection velocity 
should be of the order of Va ■ This sort of reconnection rate has never been reported to be attainable 
within plasma reconnection simulations (see [16]). However, this is the reconnection rate that is 
expected for trans-Alfvenic turbulence within the LV99 model. 


(c) Reconnection diffusion and star formation 

As we have argued earlier at length, standard flux-freezing breaks down at every point and 
time in a turbulent plasma. In that case, the only objectively meaningful way to give a 
magnetic field-line an identity over time is by tagging it with a certain plasma fluid element. 
As suggested by Axford [147], we understand the crucial feature of magnetic reconnection to 
be the "disconnection" of fluid elements that start on the same field line. The right panel of 
Figure 11 below uses data from the JHU MHD turbulence database archive to illustrate how an 
initial magnetic field line changes its connections to plasma fluid elements over time. The figure 
shows shows an initial magnetic field line, in black, decorated with eleven plasma fluid elements, 
indicated by various colors. The plasma elements are then evolved with the fluid velocity for 
about one large-eddy turnover time (t = 2.00 in imits of the simulation). The magnetic field-lines 
threading these later plasma elements are drastically different. Indeed, the plasma has "drifted" 
to distinct lines separated by distances of order the magnetic integral length (0.35 in the units 
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of the simulation). This drift occurs even though the conductivity of the simulation is high and 
the Ohmic electric fields are tiny, because their small direct effects are greatly magnified by the 
turbulence. 




Figure 11. Turbulent Splitting of a Magnetic Field Line. An initial magnetic field line in black is decorated 
with plasma elements, which are allowed to move with the fluid velocity. The magnetic field-lines that thread these plasma 
elements at the later time are drastically different, effectively “splitting" the original line. Left panel. Reconnection event 
within 3D driven turbulence. The structure of the event corresponds to the LV99 picture and is consistent with 
reconnection events studied in solar wind. 

The violation of flux freezing means that the astrophysical theories based on the concept of 
flux freezing must be revised. In particular, the standard star formation theory assumes that flux 
freezing is being violated in the partially ionized gas only due to the relative drift of neutrals and 
ions. As we discussed in sections 4(a) and 4(b) in the partially ionized gas the important magnetic 
flux violation arises from magnetic diffusion induced by turbulence. This process that was termed 
"reconnection diffusion" was identified and described in [148] (see also [106]) and successfully 
tested in the subsequent publications for the case of molecular clouds and protostellar disks, e.g. 
[68,149-152]. A comprehensive review dealing with reconnection diffusion is presented in [153]. 

The left panel illustrates magnetic the evolution of magnetic flux during the process of 
reconnection diffusion of magnetic flux out of the circumstellar accretion disk. The magnetic field 
lines are smoothed in the picture to illustrate the evolution of the mean magnetic field. 

The theory of transporting matter in turbulent magnetized medium is discussed at length 
in [154] and [153] and we refer our reader to these publications. The process was termed 
"reconnection diffusion" to stress the importance of reconnection in the the diffusive transport. 

The peculiarity of reconnection diffusion is that it requires nearly parallel magnetic field lines 
to reconnect, while the textbook description of reconnection is usually associated with anti¬ 
parallel description of magnetic field lines. One should understand that the configuration shown 
in Figure 2 is just a cross section of the magnetic fluxes depicfing the anti-parallel components of 
magnetic field. Generically, in 3D reconnection configurations the sheared component of magnetic 
field is present. The process of reconnection diffusion is closely connected with the reconnection 
between adjacent Alfvenic eddies (see upper panel of Figure 12). As a result, adjacent flux tubes 
exchange their segments with entrained plasmas and flux tubes of different eddies get connected. 
This process involves eddies of all the sizes along the cascade and ensures fast diffusion which 
has similarities with turbulent diffusion in ordinary hydrodynamic flows. The lower panel 

The lower panel of Figure 12 illustrates magnetic the evolution of magnetic flux during the 
process of reconnection diffusion of magnetic flux out of the circumstellar accretion disk. The 
magnetic field lines are smoothed in the picture to illustrate the evolution of fhe mean magnetic 


field. 


Reconnection diffusion should nof necessarily be understood as a concept that makes the 
earlier theories of sfar formation invalid. In fact, turbulence in dark cores giving birth to 
stars may be reduced and this may make the traditional ambipolar diffusion, i.e. the drift of 
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Figure 12. Upper panel: Reconnection diffusion: exchange of flux with entrained matter. Illustration of the mixing of 
matter and magnetic fields due to reconnection as two flux tubes of different eddies interact. Only one scale of turbulent 
motions is shown. In real turbulent cascade such interactions proceed at every scale of turbulent motions. From [154]. 
Lover panel: Reconnection diffusion in in an accretion disk. Illustration of the process using smoothen lines. 
From Casanova et al. (2015). 


neutrals in relation to ions, important. However, reconnection diffusion must be a part of star 
formation paradigm. In fact, it can successfully explain many pieces of observational data that 
are completely puzzling within the ambipolar diffusion paradigm. This includes, for instance, the 
famous "magnetic breaking catastrophe" for accretion disks which is the inability of removing 
magnetic flux from accretion disks fast enough to enable forming of such disks. Similarly, the 
poor correlation of density and magnetic field in interstellar media is also impossible to explain 
on the basis of ambiplar diffusion. A comprehensive review dealing with reconnection diffusion 
is presented in [153]. Closely related is the recent development of the "turbulent general magnetic 
reconnection" (TGMR) theory in [35]. The starting point of this theory is the understanding that 
magnetic field-line "motion" can be objectively defined only relative to plasma fluid elements 
and their magnetic connectivity. In star formation for example, the magnetic field-lines threading 
the protostar will appear to "slip" relative to the ambient ISM, whereas the field-lines embedded 
in the ISM will likewise appear to "slip" through the collapsing magnetic cloud. Neither picture 
is more correct than the other and, indeed, one cannot uniquely define a "motion" of the field¬ 
lines. However, it has been shown in [35] that the lines wandering between the protostar and 
the surrounding ISM acquire a imique slip-velocity per unit arc-length of field-line, which is 
completely independent of which end of the line is regarded to be the "foot-point" tied to the 
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plasma. This so-called slip-velocity source is given by the expression 
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E = -{\7 xR)^/\B\ (7.6) 

where R is the non-ideal electric field in the generalized Ohm's law E + u x B = R for the plasma, 
and "_L" denotes the component perpendicular to B. It is only by spatially wandering/intersecting 
a region with non-vanishing E that a field-line can evade the "frozen-in" condition. Furthermore, 
in a turbulent plasma the slip-source E is enormous, even though the electric fields R are tiny 
and the plasma nearly ideal. This approach is essentially a refinement of the LV99 idea that 
field-lines must wander into microscopic "current sheets" in order to break the flux-freezing 
constraint. [35] applied this theory to explain observed deviations from the Parker spiral model 
of the interplanetary magnetic field in our own solar system, due to "slippage" of the field-lines 
through the turbulent solar wind. 

(d) Solar flares and gamma ray bursts 

The picture of flares of reconnection described in section 4d is broadly supported by current 
observations and numerical simulations of solar flares and CME's. For example, simulations 
by [155] of the "breakout model" of CME initiation show that an extremely complex magnetic 
line structure develops in the ejecta during and after the initial breakout reconnection phase, 
even under the severe numerical resolution constraints of such simulations. In the very high 
Lundquist-number solar environment, this complex field must correspond to a strongly turbulent 
state, within which the subsequent "anti-breakout reconnection" and post-CME current sheet 
occur. Direct observations of such current sheets [124,156] verify the presence of strong turbulence 
and greatly thickened reconnection zones, consistent with the LV99 model. In the numerical 
simulations, the "trigger" of the initial breakout reconnection is numerical resistivity and there 
is no evidence of turbulence or complex field-structure during the eruptive flare onset. This is 
very likely to be a result of the limitations on resolution, however, and we expect that developing 
turbulence will accelerate reconnection in this phase of the flare as well. 

While the details of the physical processes discussed above can be altered, it is clear that LV99 
reconnection induces bursts in highly magnetized plasmas. This can be applicable not only to the 
solar environment but also to more exotic environments, e.g. to gamma ray bursts. The model 
of gamma ray bursts based on LV99 reconnection was suggested in [115]. It was elaborated and 
compared with observations in [36]. Currently, the latter model is considered promising and it 
attracts a lot of attention of researchers. Elares of reconnection that we described above can also 
be important for compact sources, like pulsars and black holes in microquasars and AGNs [116]. 
We would like to note that LV99 reconnection is getting more applications related to emission of 
astrophysical objects. Eor instance, recently it has been discussed to explain the radio and gamma 
ray emission arising through accretion on black holes [157] as well as for describing the radiation 
of microquasars [158]. 

(e) Turbulent reconnection and particle acceleration 

Turbulent reconnection provides the way of the first order Eermi acceleration as it is illustrated 
in Eigure 13. The efficiency of the process is ensured by LV99 model being the volume-filling 
reconnection®. 

The left panel of Eigure 13 illustrates a situation when the particle anisotropy which arises 
from particles preferentially accelerated in direction parallel to magnetic field. Similarly, [159] 
showed fhat the first order Fermi acceleration can also happen in terms of the perpendicular to the 
magnetic field component of particle momentum. This is illustrated in the right panel of Figure 13. 
There fhe particle with a large Larmour radius is bouncing back and forth between converging 

®We would like to stress that Figure 2 exemplifies only the first moment of reconnection when the fluxes are just brought 
together. As the reconnection develops the volume of thickness A becomes filled with the reconnected 3D flux ropes moving 
in the opposite directions. 
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Figure 13. Left Panel: Cosmic rays spiral about a reconnected magnetic field line and bounce back at points A and B. 
The reconnected regions move towards each other with the reconnection velocity Vr. From [148]. Right Panel: Particles 
with a large Larmor radius gyrate about the magnetic field shared by two reconnecting fluxes (the latter is frequently 
referred to as “guide field”. As the particle interacts with converging magnetized flow corresponding to the reconnecting 
components of magnetic field, the particle gets energy gain during every gyration. From [159]. 


mirrors of reconnecting magnetic field systematically getting an increase of the perpendicular 
component of its momentum. Both processes take place in reconnection layers. 

Disregarding the backreaction one can get the spectrum of accelerated cosmic rays [116,148]: 

N{E)dE = constiE~^/‘^dE, (7.7) 

This result is the result of acceleration in the absence of compression (see [160]). The first order 
acceleration of particles entrained on the contracting magnetic loop can be understood from the 
Liouville theorem. In the process of fhe magnetic tubes contraction a regular increase of the 
particle's energies is expected. The requirement for the process to proceed efficiently is to keep 
the accelerated particles within the contracting magnetic loop. This introduces limitations on the 
particle diffusivity perpendicular to the magnetic field direction. The process in Figure 13 (left 
panel) was discussed in [11] in relation to the acceleration of particles in collisionless reconnection. 
There by accounting for the backreaction of particles the authors obtained a more shallow 
spectrum. 

Testing of particle acceleration in turbulent reconnection was performed in [34] and its results 
are presented in Figure 14. The Figure 14 shows the evolution of fhe kinetic energy of the particles. 
After injection, a large fraction of test particles accelerates and the particle energy growth occurs 
(see also the energy spectrum at t = 5 in the detail at the bottom right). This is explained by a 
combination of two effects: the presence of a large number of converging small scale current sheets 
and the broadening of the acceleration region due to the turbulence. The acceleration process is 
clearly a first order Fermi process, and involves larger number of particles, since the size of the 
acceleration zone and the number of scatterers naturally increases by the presence of turbulence. 
Moreover, the reconnection speed, which in this case is independent of resistivity [31,33] and 
determines the velocity at which the current sheets scatter particles, has been naturally increased 
as well (i.e. Vrec ~ Va). During this stage the acceleration rate is in the range 2.48 — 2.75. 

The process of acceleration via turbulent reconnection is expected to be widespread. In 
particular, it has been discussed in [135] as a cause of fhe anomalous cosmic rays observed by 
Voyagers and in [161] as a source of the observed cosmic ray anisotropies. We expect turbulent 
reconnection to accelerate energetic particles in relativistic environments, like those related to 
accretion disks and relativistic jets as well as in gamma ray bursts (see [21] for a review). The 
latter process discuss first in [115], has been given strong observational support in [36]. 
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Time 

Figure 14. Particle kinetic energy distributions for 10,000 protons injected in the fast magnetic reconnection domain. The 
colors indicate which velocity component is accelerated (red or blue for parallel or perpendicular, respectively). The energy 
is normalized by the rest proton mass. Subplot shows the particle energy distributions at t = 5.0. Models with Bqz = 0.1, 
rj = 10“®, and the resolution 256x512x256 is shown. 


8. Connection with other ideas of fast reconnection 

As we mentioned in the Introduction the reconnection research is a vast vigorously developing 
field and is not limited to reconnection in MHD approximation in turbulent fluids that we deal 
in this review. A lot of experimental research is done for the Earth magnetosphere and laboratory 
plasmas, where the MHD description is not valid. Some plasmas may not be turbulent either. 
Below we briefly outline the connection of our turbulent model with other directions of research. 

(a) Plasmoid/tearing mode reconnection 

Plasmoid/tearing mode recormection is currently a vibrant direction of reconnection research 
(see [23,162]). The work in this direction shows that SP reconnection is unstable for sufficiently 
large Ludquist numbers. What should be kept in mind that in 3D the thicker outflows induced 
by plasmoid/tearing reconnection inevitably induce turbulence. This corresponds to both PIC 
and MHD simulations that we discussed earlier. This also inevitable on the theoretical grounds. 
Indeed, from the mass conservation constraint requirement in order to have fast reconnection one 
has to increase the outflow region thickness in proportion to Lx, which means the proportionality 
to the Lundquist number S. The Reynolds number Re of the outflow is AVa/v, where v is 
viscosity, grows also as S. The outflow gets turbulent for sufficiently large Re. It is natural to 
assume that once the shearing rate introduced by eddies is larger than the rate of the tearing 
instability growth, the instability should get suppressed. 

If one assumes that tearing is the necessary requirement for fast recormection this entails the 
conclusion that tearing should proceed at the critically damped rate, which implies that the Re 
number and therefore A should not increase. This entails, however, the decrease of reconnection 
rate driven by tearing in proportion ~ 5 as a result of mass conservation. As a result, the 
reconnection should stop being fast. Fortunately, we know that turbulence itself provides fast 
reconnection irrespectively whether tearing is involved or not. Thus one may conclude that the 
tearing reconnection, similar to the SP recormection, should be applicable to a limited range of 
S for realistic magnetized plasmas with low viscosity in the perpendicular to magnetic field 
direction. 

Tearing may be important for initiating turbulence and transiting from the laminar initial state. 
To what extend tearing is required is not clear from the 3D simulations that we discussed above. 
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Those suggest the importance of Kelvin-Helmholdz instability, but whether tearing plays any 
role must still be explored. However, if reconnection was excited by tearing/plasmoid instability 
generically we expect a transition to the regime of turbulent reconnection. 

Another limitation to the applicability of tearing reconnection arises from its speed. The 
reported rates do not exceed a small fraction of Alfven speed. However, as we discuss later 
magnetic turbulence requires reconnection speeds which are substantially larger than that value 
for the theory of MHD turbulence to be self-consistent (LV99, ELVll). In this situation, the 
dominance of turbulent reconnection seems inevitable. 

(b) Reconnection in 2D in presence of turbulence 

[42,163] explored numerically turbulent reconnection in 2D. As a theoretical motivation the 
authors emphasized analogies between the magnetic reconnection layer at high Limdquist 
numbers and homogeneous MHD turbulence. They also pointed out various turbulence 
mechanisms that would enhance reconnection rates, including multiple X-points as reconnection 
sites, compressibility effects, motional EMF of magnetic bubbles adverting out of the reconnection 
zone. However, the authors did not understand the importance of "spontaneous stochasticity" of 
field lines and of Lagrangian trajectories and they did not arrive at an analytical prediction for 
the reconnection speed. Although an enhancement of the reconnection rate was reported in their 
numerical study, but the setup precluded the calculation of a long-term average reconnection rate. 

The relation of this study with LV99 is not clear, as the nature of turbulence in 2D is different. 
In particular, shear-Alfven waves that play the dominant role in 3D MHD turbulence according 
to GS95 are entirely lacking in 2D, where only pseudo-Alfven wave modes exist. We believe that 
the question whether turbulence is fast has not been resolved yet if we judge from the available 
publications. For instance, in a more recent study along the lines of the approach in [163], i.e. 
in [164], the effects of small-scale turbulence on 2D reconnection were studied and no significant 
effects of turbulence on reconnection was reported. ]165] have more recently made a study of 
Ohmic electric fields at X-points in homogeneous, decaying 2D MHD turbulence. However, they 
studied a case of small-scale magnetic reconnection and their results are not directly relevant to 
the issue of reconnection of large-scale flux tubes that we deal in this review. 

Instead of studying bulk reconnection in 2D turbulence as the aforementioned studies did, 
[166] and [167] studied large scale reconnection^, which is advantageous if the determination 
of fhe actual recormection rates is sought. The two groups reached different conclusions On the 
one hand, [166] had a better resolution but used periodic boundary conditions, which strongly 
constrain the ability to do averaging of the reconnection rate and the attainment of the steady state 
for reconnection. They inferred from their data that the 2D turbulent reconnection rate may be 
independent of resistivity. On the other hand, [167] used smaller data cubes but longer averaging, 
which is enabled by their outflow boundary conditions. They concluded that the reconnection 
does depend on resistivity and therefore is slow. 

In view of the difference of MHD turbulence in 2D and 3D we do not view the recormection 
studies in 2D turbulence as directly relevant in any astrophysical settings. Even if eventually 2D 
reconnection is proven to be fast, the reconnection rate is expected to have different dependences 
on turbulent power. 


(c) Turbulent reconnection models based on mean field approach 

Guo et al. ]44] modified and extended ideas originally proposed in [102] and suggested their 
model of fast turbulent reconnection. Both papers use mean field approach, but unlike the study 
in [102] which concluded that turbulence cannot accelerate reconnection the more recent study 
obtains expressions for fast reconnection. These expressions are different from those in LV99 and 

^The enhancement of 2D large scale reconnection was reported starting from 2007 at a few conferences by the authors of [33], 
but for them the 2D study was a testing ground for the realistic 3D simulations to test LV99. Thus these results were never 
published. 
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seem to grossly contradict to the numerical testing of turbulent reconnection in [33]. Another 
model of turbulent reconnection based on the mean field approach is presented in [168] and it 
also suffers with the problems of using the mean field approach for reconnection that we describe 
below. 

The mean field approach invoked in the aforementioned studies is plagued by poor 
foundations and conceptual inconsistencies, however [32]. In such an approach effects of 
turbulence are described using parameters such as anisotropic turbulent magnetic diffusivity 
experienced by the fields once averaged over ensembles. The problem is that it is the lines of the 
full magnetic field that must be rapidly reconnected, not just the lines of the mean field. ELVll 
sfress that the former implies the latter, but not conversely. No mean-field approach can claim 
to have explained the observed rapid pace of magnetic reconnection unless it is shown that the 
recormection rates obtained in the theory are strictly independent of the length and timescales of 
the averaging. Naturally, it is impossible to get reliable results applying mean field approach to 
recormection (see more discussion in ELVll). 

Other attempts to get fast magnetic recormection from turbulence are related to the so-called 
hyper-resistivity concept [41,43,169,170], which is another attempt to derive fast recormection 
from turbulence within the context of mean-field resistive MHD. Apart from the generic problems 
of using the mean field approach, we would like to point out that the derivation of the hyper¬ 
resistivity is questionable from a different point of view. The form of the parallel electric field 
is derived from magnetic helicity conservation. Integrating by parts one obtains a term which 
looks like an effective resistivity proportional to the magnetic helicity current. There are several 
assumptions implicit in this derivation, however. Fundamental to the hyper-resistive approach is 
the assumption that the magnetic helicity of mean fields and of small scale, statistically stationary 
turbulent fields are separately conserved, up to tiny resistivity effects. However, this ignores 
magnetic helicity fluxes through open boundaries, essential for stationary recormection, that 
vitiate the conservation constraint. 

As we discuss further, a common misunderstanding is that "resistivity arising from 
turbulence" is a real plasma non-ideality "created" by the turbulence. However, such apparent 
non-ideality is strongly dependent on the length and timescales of the averaging. It appears only 
as a consequence of observing the plasma dynamics at a low resolution, so that the coarse-grained 
velocity and magnetic field that are observed will no longer satisfy fhe microscopic equations of 
motion. This coarse-graining or averaging is a purely passive operation which doesn't change the 
actual plasma dynamics but only corresponds to "taking off one's specfacles". It is clear that one 
cannot create true, physical non-ideal electric fields by removing one's eyeglasses! Such apparent 
non-ideality in a turbulent plasma observed at length-scales in the inertial-range or larger is a 
valid representation of the effects of turbulent eddies at smaller scales. However, such apparent 
non-ideality is not accurately represented by an effective "resistivity", a representation which in 
the fluid turbulence literature has been labelled the "gradient-transport fallacy" [171]. It is also 
clear that no mean-field or coarse-graining approach can claim to have explained the observed 
rapid pace of magnetic reconnection imless it is shown that the recormection rates obtained in the 
theory are strictly independent of fhe length and timescales of the averaging [35,172]. 

More detailed discussion of fhe conceptual problems of the hyper-resistivity concept and mean 
field approach to magnetic reconnection is presented in [95] and ELVll. 

(d) Indirect evidence for turbulent reconnection 

A study of tearing instability of current sheets in the presence of background 2D turbulence that 
observed the formation of large-scale islands was performed in [173]. While one can argue that 
observed long-lived islands are the artifact of adopted 2D geometry, the authors present evidence 
for fast energy dissipation in 2D MHD turbulence and show that this result does not change as 
they change the resolution. More recently [174] provided evidence for fast dissipation also in 3D 
MHD turbulence. This phenomenon is consistent with the idea of fast reconnection, but cannot be 
treated as a direct evidence of fhe process. Indeed, fast dissipation and fast magnetic recormection 
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are rather different physical processes, dealing with decrease of energy on the one hand and 
decrease of magnefic flux on the other. 

Works by Galsgaard and Nordlund, e.g. [62], could also be interpreted as an indirect support 
for fast reconnection. The authors showed that in their simulations they could not produce highly 
twisted magnetic fields. One possible inferpretation of this result could be the fast relaxation 
of magnetic field via reconnecfion^*^ However, in view of many uncerfainties of fhe numerical 
studies, this relation is unclear. With low resolution involved in the simulations the Reynolds 
numbers could not allow a turbulent inertial-range. 

9. Concluding Remarks 

(a) Turbulent reconnection and “turbulent resistivity” 

As we discussed in the review, the violation of flux freezing and diffusivity of magnetic field 
that contradicts to the Alfven theorem follows from the LV99 model of fast reconnection. This, 
however, is sometimes, misunderstood as our using some sort of 'Turbulent resistivity". As we 
mentioned in the review, this confusion is common for many papers. Therefore we discuss this 
issue here in more detail. It is possible to show that "turbulent resistivity" description has fatal 
problems of inaccuracy and unreliabilify, due to its poor physical foundations for turbulent flow. 
It is true that coarse-graining the MHD equations by eliminating modes at scales smaller than 
some length I will introduce a "turbulent electric field", i.e. an effective field acting on the large 
scales induced by motions of magnetized eddies at smaller scales. However, it is well-known 
in the fluid dynamics community that the resulting turbulent transport is not "down-gradient" 
and not well-represented by an enhanced diffusivity. The physical reason is that turbulence lacks 
the separation in scales to justify a simple "eddy-resistivity" description. As a consequence, 
energy is often not absorbed by the smaller eddies, but supplied by them, a phenomenon called 
"backscatter". In magnetic reconnection, the turbulent electric field offen creates magnetic flux 
rafher than destroys it. 

If we know the recormection rate, e.g. from LV99, then an eddy-resistivity can always be tuned 
by hand to achieve that rate. But this is engineering, not science. While the tuned reconnection 
rate will be correct by construction, other predictions will be wrong. The required large eddy- 
resistivity will smooth out all turbulence magnetic structure below the coarse-graining scale 
L In reality, the turbulence will produce strong small-scale inhomogeneities, such as current 
sheets, from the scale I down to the micro-scale. In addition, field-lines in the flow smoothed 
by eddy-resistivity will not show the explosive, super-diffusive Richardson-type separation at 
scales below 1. These are just example of effects that will be lost if the wrong concept of 
"eddy resistivity" is adopted. Note, that the aforementioned are important for understanding 
particle transport/scattering/acceleration in the turbulent recormection zone. Continuing with 
the list, we can point out that in the case of relativistic recormection, turbulent resistivities will 
introduce acausal, faster than light propagation effects. Nevertheless, the worst feature of the 
crude "eddy-resistivity" parameterization is its unreliability: because it has no sound scientific 
basis whatsoever, it cannot be applied with any confidence to astrophysical problems. Therefore 
it is pointless to talk about "turbulent resistivity" for fhe problems that we discussed in the review, 
e.g. solar flares, sfar formafion, gamma ray bursfs. 

Equivalenfly, the stochastic flux freezing [32] closely related to the fast turbulent recormection 
concept is definitely not equivalent to the dissipation of magnetic field by resistivity . While the 
parametrization of some particular effects of turbulent fluid may be achieved in models with 
different physics, e.g. of fluids wifh enormously enhanced resisfivify, the difference in physics 
will inevifably result in other effects being wrongly represented by this effect. For instance, 
turbulence with fluid having resistivity corresponding to the value of "turbulent resistivity" 
must have magnetic field and fluid decoupled on most of its inertia range turbulent scale, 

this case, these observations could be related to the numerical finding of [175] which shows that reconnecting magnetic 
configurations spontaneously get chaotic and dissipate, which, as discussed in [176], may be related to the LV99 model. 
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i.e. the turbulence should not be affected by magnetic field in gross contradiction with theory, 
observations and numerical simulations. Magnetic helicity conservation which is essential for 
astrophysical dynamo should also be grossly violated^ 

The approach advocafed by us in discussing furbulenf reconnection is quite different. It is not 
based on coarse-graining. The spontaneous stochasticity of magnetic field-lines and of Lagrangian 
frajecfories (plasma fluid element histories) is a real, verified physical phenomenon in turbulent 
fluids. Whereas "eddy-resistivity" ideas predict that magnetic flux is destroyed by turbulence, 
our work shows that turbulent spontaneous stochasticity transforms magnetic-flux conservation 
into a stochastic conservation law. Because spontaneously stochastic world-lines in relativistic 
turbulence must remain within the light-cone, no acausal effects such as produced by "eddy- 
resistivity" will be predicted. Our approach is based on fundamental scientific progress in the 
understanding of furbulence, not on engineering parameterizations. 

(b) Goldreich-Sridhar turbulence and turbulent reconnection 

GS95 turbulence is a theory accepted by a substantial part of the astrophysical community 
(see [89,91] for reviews). Born oufside the mainstream community of turbulence experts, it was 
initially mostly ignored^^ but then was accepted under the pressure of numerical results. As we 
mentioned, the debates are not settled about the validity of possible modifications of fhe model. In 
parallel, some part of the community is still using the so-called 2D plus slab model of furbulence 
(see [177]) in spife of the fact that it has no support via numerical simulations with isotropic 
driving. We consider the latter as some parametrization of fhe actual heliospheric turbulence over 
a limited range. This parametrization is not physically or numerically motivated and therefore is 
not considered within the reconnection model. 

The LV99 and our subsequent studies mentioned in the review employed GS95 model. 
However, we would like to stress again that none of our principal resulfs on fast turbulent 
reconnection and the physics of turbulent reconnection will be changed if instead of GS95 any 
other existing model of strong MHD turbulence is used, provided that this model satisfies the 
constraints that are given by the existing numerical simulations. The corresponding expressions 
for reconnection rafes obtained a wide variety of turbulence indexes are provided in LV99. 
At the same time, we discussed in the review that fast LV99 reconnection makes GS95 model 
self-consistent. 

(c) 2D and 3D reconnection 

Numerical simulations are very demanding in 3D and therefore as therefore the numerical 
research attempts initially to attack the problem of reduced dimensions. In terms of recormection 
atfempts to attack the problem with 2D simulations are widely spread. While 2D simulations can 
get insight to some processes, the relation between 2D and 3D reconnection is far from trivial. For 
instance, from general theoretical positions the importance of 3D for reconnection was advocated 
by [178,179]. 

In general, a radical change of physics related with the use of the 2D instead of actual 3D 
is very common for complex physical systems and the problem of obtaining misleading results 
extrapolating those obtained in 2D for the actual 3D systems goes beyond turbulence. Every time 
when the 2D physics is employed, it is essential to prove that the results stay the same in 3D. This, 
for instance, has not been done in the case of 2D turbulent reconnection [42] and we believe that 
this may not be possible due to fundamental differences of turbulence physics (see ELVll). 

Even in the case when 2D reconnection reflects the physics common to 3D, e.g. in the case of 
fearing reconnection (e.g. [23]), we claim that the development of the instability in 3D and 2D 

^^This is a serious mistake of a number of numerical simulations of galactic dynamos where to "simulate" effects of turbulent 
diffusion the Ohmic resistivity y of the order of Vl L is used. Surely these simulations do not represent the actual fast 
astrophysical dynamo, but only a slow one. 

^^The enthusiasm of accepting alternative theories that, e.g. provide more traditional for the MHD community Kraichnan 
index of —3/2 [81,82], may also be partially explained by this fact. 


rspa.royalsocietypublishing.org Proc R Soc A 0000000 




may be very different. The 3D configurations are more prone to secondary instabilities and to the 
development of the fully turbulent state in which the initial instabilities may not be dominant or 
even important. 


36 


(d) Turbulent reconnection and plasma effects 

A substantial part of the reconnection research is based on exploring plasma physics effects on 
reconnection (see [16,180] for reviews). LV99 shows that reconnection rates should not depend on 
plasma microphysics in the presence of turbulence. This conclusion was supported by numerical 
study in [33] where plasma effects were simulated by introducing anomalous resistivity. The 
subdominance of reconnection arising from the Hall effect to that arising due to turbulence 
was shown analytically in ELVll. A more rigorous comparison of the turbulence induced 
recormection with that induced by other terms in the Generalized Ohm Equation was provided in 
[35] where it was shown that for typical astrophysical parameters turbulence effects are absolutely 
dominant. 

Nevertheless, the studies that show that magnetic reconnection can be fast in the absence of 
turbulence (see [16] for a review and ref. therein) poses interesting questions on the actual role 
of turbulence. There are, for instance, suggestions that tearing of the current sheet may make 
collisionless plasma effects applicable to magnetic reconnection on large astrophysical scales (see 
[29] for a review). Clearly in plasma one should consider both small scale plasma turbulence as 
well as large scale turbulence that obeys MHD treatment. In addition, the very issue of plasma 
collisionality is frequently unclear. Indeed, apart from Coulomb collisions, ions may be scattered 
by magnetic inhomogeneities that arise due to a number of instabilities (e.g. firehose, mirror) in 
collisionaless plasmas and this may make plasma effectively collisional (see [51,69]) in agreement 
with recent simulations ]68]. In this situation, the MHD description should be applicable even to 
plasmas which are formally collisionless. 

Our arguments in the review suggest that plasma effects cannot dominate large scale 
astrophysical reconnection in the presence of turbulence. However, we accept that the issue is 
a subject of interesting debates and more testing is valuable. 


(e) Present state of turbulent reconnection theory and outstanding 
questions 

We would like to emphasize that at present the turbulent reconnection theory does not amount 
to the LV99 model only. It is also ELVll where the LV99 expressions were reproduced using a 
very different approach that follows from the recent advances of the Lagrangian description of 
turbulence. It is also a very recent paper by [35] where the effects of turbulence were included 
within the Generalized Ohm's law and were shown to be in agreement with results of the two 
approaches above. The theoretical foundation of turbulent reconnection got substantial support 
from the recently developed concept of "spontaneous stochasticity" [32]. 

The predictions of the turbulent reconnection have been successfully tested both with direct 
simulations of turbulent reconnection layer [33,34] and through the violation of flux conservation 
that the turbulent reconnection entailed [122]. As we discuss in the review, more promising 
numerical tests of reconnection, including turbulence being self-driven are under way. 

The turbulent reconnection has shown promise in explaining various astrophysical problems, 
as well as in addressing problems facing solar physics and heliospheric research. Some of these 
are discussed in this review, while others are discussed in specialized reviews [21,153,181]. 

At the same time, a number of questions remains not answered. It is obvious that LV99 model 
is a very simplified model. It does not take into account many effects, e.g. effects of plasma 
compressibility, turbulence intermittency, velocity and magnetic field shear. To obtain analytical 
results it assumes the turbulence is presented by a single power law and disregards the deviations 
arising from multiple scales of energy injection. Ohmic and viscous dissipation etc. It deals with 


rspa.royalsocietypublishing.org Proc R Soc A 0000000 




37 


isothermal MHD description of the process and does not account for relativistic effects. The role 
of collisionless plasma effecfs in turbulent reconnection is hotly debated. 

These limitations of the model are gradually dealt with. For instance, a modification of our 
understanding of GS95 cascade was discussed in our review in relation to describing magnetically 
dominated perturbations arising from magnetic reconnection. We also discussed accounting for 
the partial ionization of plasmas. Role of plasma effects is also being clarified (see [35] and ref. 
therein). Nevertheless, we are at the very beginning of our sfudies of turbulent reconnection 
and its consequences. Therefore we expect many surprises and discoveries on the way to fully 
understanding of fhe intricate relation of turbulence and reconnection. 

(f) Turbulence as a converging point for reconnection research 

LV99 model is the one that describes the dynamics of reconnection within turbulent fluids in 
MHD regime. By itself, the study of dynamics of magnetic fields in MHD regime, irrespectively to 
any plasma physics is a well motivated direction. However, astrophysical environments are filled 
with turbulent plasmas. Thus for astrophysical applications it is important to define the domain of 
applicability of turbulence-based versus plasma-based reconnection. First of all, we should stress 
that there is no single mode of reconnection. Astrophysical and laboratory environments present 
an extensive variety of conditions for magnetic reconnection. 

Numerical experiments show that laminar Sweet-Parker reconnection is feasible in the regime 
of low Ludquist numbers, then the laminar picture fails and tearing gets important. As we further 
increase the length of the reconnection sheet, even without external driving, in 3D low viscosity 
plasma the transition to turbulence is inevitable. Whether this transfers fhe reconnection to purely 
turbulent reconnection or plasma effects are important for determining the reconnection rate 
when the outflow is fully turbulent is a subject of ongoing debates. In our work we provided 
arguments in support to the former solution, i.e. that the transition to the turbulent state when 
plasma effect do not change the reconnection rate is most relevant for mosf asfrophysical settings. 
This does not exclude that in some particular circumstances, e.g. for the onset of turbulent 
reconnection from initially laminar state, for current sheet which thickness is comparable with 
ion inertial length, as this is the case of magnetoshpere^^, the plasma effects are important. For 
other opinions and outgoing debates we refer our reader to [16]. Below we, however, point out to 
the tendency that the reconnection research has demonstrated. 

Recent years have demonstrated the convergence of turbulent reconnection in LV99 and 
other directions of reconnection research. For instance, models of collisionless reconnection have 
acquired several features in common with the LV99 model. In particular, they have moved to 
consideration of volume-filling reconnection (see [11]). While much of fhe discussion may still 
be centered around 2D magnetic islands produced by reconnection, in three dimensions these 
islands are expected to evolve into contracting 3D loops or ropes [182] introducing stochasticity 
to the reconnection zone. Moreover, it is more and more realized that the 3D geometry of 
reconnection is essential and that the 2D studies may be misleading. 

The departure from the concept of laminar reconnection and the introduction of magnetic 
stochasticity is also apparent in a number of recent papers appealing to the tearing mode 
instability to drive fast reconnection (see [23,27]). These studies showed that tearing modes 
do not require collisionless environments and thus collisionality is not a necessary ingredient 
of fast reconnection. Finally, the development of turbulence in 3D numerical simulations of 
reconnection (see section 4c) clearly testifies that the reconnection induces turbulence even if the 
initial reconnection conditions are laminar. 

All in all, in the last decade, the models competing with LV99 have undergone a substantial 
evolution, from 2D collisionless reconnection based mostly on Hall effect to 3D reconnection 
where the collisionless condition is no more required, Hall effect is not employed, but 
magnetic stochasticity and turbulence play an important role in the thick reconnection regions. 

^^Plasma turbulence may be still important for such reconnection, but this type of turbulent recormection is not described by 
LV99 model. 
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Nevertheless, we want to stress that collisionless reconnection may be suitable for the description 
of reconnection when the reconnecting flux-structures are comparable with the ion gyro scale, 
which is the case of the reconnection studied situ in the magnetosphere. However, this is a special 
case of magnetic reconnection which makes its very atypical generic astrophysical settups where 
reconnection involves scales many orders of magnetude of the gyroradius involved. Even in this 
case we may expect the development of turbulence, but this would not be MHD turbulence 
which makes LV99 theory not applicable to it. Conversely, it can be shown by exact analytical 
estimates [35] that the direct effect of fhe microscopic plasma non-idealities are negligible for 
reconnection at scales vastly larger than ion gyroradius. 
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